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Kaon polarization operator in dense baryonic matter of arbitrary isotopic composition is calculated 
including s- and p-wave kaon-baryon interactions. The regular part of the polarization operator 
is extracted from the realistic kaon-nucleon interaction based on the chiral and 1/Nc expansion. 
Contributions of the A(1116), E(1195), E*(1385) resonances are taken explicitly into account in the 
pole and regular terms with inclusion of mean-field potentials. The baryon-baryon correlations are 
incorporated and fluctuation contributions are estimated. Results are applied for K~ in neutron 
star matter. Within our model a second-order phase transition to the s-wave K~ condensate state 
occurs at pc ~ 4po once the baryon-baryon correlations are included. We show that the second- 
order phase transition to the p-wave K~ condensate state may occur at densities Pc ~ 3 5po 
in dependence on the parameter choice. We demonstrate that a first-order phase transition to a 
proton-enriched (approximately isospin-symmetric) nucleon matter with a p-wave K~ condensate 
can occur at smaller densities, p 2po. The transition is accompanied by the suppression of hyperon 
concentrations. 
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I. INTRODUCTION 

Strangeness modes in compressed hadronic matter stay 
in the focus of interest over the last decade. Strangeness 
is considered as a good probe for dynamics of heavy-ion 
collisions. The vast amount of data is accumulated for 
different collision energy regimes at GANIL, GSI, GERN 
and BNL facilities [1]. New data are in advent [2]. Under¬ 
standing of strangeness production requires a systematic 
study of the evolution of virtual strangeness modes in the 
quark-gluon plasma and in the soup of virtual hadrons. 
At break-up stage these virtual modes are redistributed 
between real strange particles which can be observed in 
experiment. 

Other interesting topic is a strangeness content of neu¬ 
tron stars. With density increase strangeness can be lib¬ 
erated and face up in the hlling of the hyperon Fermi seas 
and/or in the creation of kaon condensates. Both men¬ 
tioned topics need a better knowledge of the kaon-baryon 
interaction in dense baryonic matter. Question on princi¬ 
pal possibility of the kaon condensation in dense nucleon 
matter was risen in Refs. [3] and [4]. Occurrence of the 
kaon condensation in neutron star interiors may have in¬ 
teresting observational consequences: (i) The softening 
of the equation of state (EoS), due to the appearance of 
a kaon condensate phase, lowers the maximum neutron 
star mass and could induce transformation of neutron 
stars into low-mass black holes [5]. (ii) Kaon conden¬ 
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sation is predicted to be accompanied with the change 
of the nucleon isospin composition from the neutron- 
enriched star {N >> Z) to the ’nuclear’ star [N ~ Z) 
or even to the ’proton’ star {N < Z), where the elec¬ 
tric charge of protons is compensated by the charge of 
the condensed kaons [6]. (iii) The enhanced neutrino- 
emission processes, occurring on protons in the proton- 
enriched matter and on the kaon condensate field, lead 
to substantially more rapid cooling of the star [7]. 

The K~ condensate is created in neutron stars due to 
weak multi-particle processes 

e-+X^K-+X', n + X->p + K-+X', (1) 

in which electrons are replaced by K~ mesons, and 
neutrons are converted into protons and K~ [8]. The 
symbolic writing (1) assumes that surrounding baryons 
{X, X') assure the momentum conservation, therefore the 
critical point is determined by the energy balance only. 
These processes become possible, if the electron chemical 
potential pe exceeds the minimal K~ energy, 

a;““(A;^) = min{w™(fc)} , 

k 

where a;“‘°(|fc|) is the K~ energy at the lowest (index 
” min”) quasiparticle spectrum branch of K~ excitations 
in neutron star matter. The works [4, 8, 9] and subse¬ 
quent ones have postulated that there is only one kaon 
branch, for which u;{k = 0) rriK {mK is the kaon 
mass) for the baryon density p ^ 0, and the minimum 
is achieved at km = 0. Then the critical point of the s- 
wave K~ condensation in a second-order phase transition 
is determined by the condition a;(fe = 0) = pte- The hrst- 
order phase transition to the kaon condensate state was 
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investigated in Ref. [10] applying the Maxwell construc¬ 
tion principle and in Ref. [11] according to Gibbs criteria. 
The p-wave kaon-nucleon interaction, which changes the 
momentum dependence of the kaon spectrum, was dis¬ 
regarded in those works. The p-wave A(1116)-nucleon 
hole and S(1195)-nucleon hole contributions to the kaon 
polarization operator were introduced in Ref. [12] in the 
framework of the chiral SU(3) symmetry. However then 
authors focused on the discussion of the s-wave kaon con¬ 
densation and considered the polarization operator at 
zero momentum. 

Ref. [6] worked out a possibility of the p-wave kaon 
condensation. The kaon polarization operator was con¬ 
structed with inclusion of the A(1116)-nucleon-hole and 
S(1195)-nucleon-hole contributions in the p-wave part 
of the kaon polarization operator and the kaon-pion and 
kaon-kaon interactions. The multi-branch spectrum of 
K~ mesons was found and the possibility of the p-wave 
kaon condensation related to the population of hyperon- 
nucleon-hole modes was demonstrated. Possibilities of 
first-order phase transitions in neutron star interiors to 
a proton-enriched matter with a p-wave K~ condensate 
and to a neutron-enriched matter with a p-wave con¬ 
densate were suggested. 

Ref. [13] considered the case of a large hyperon admix¬ 
ture in the neutron star core. In such a medium both K~ 
and K'^ spectra possess extra branches associated with 
the particle-hole excitations S — A“^, S — for K~ 
and N — A”^, N — A“^ for Ar+ (hole states are labeled 
here by tense (-1)). At quite large kaon momenta the 
branches of and K~ spectra merge, that signals an 
instability with respect to K~ pair creation. 

The analyses [6, 13] relied heavily on the pole approx¬ 
imation for the particle-hole diagrams. The final width 
effects were thereby disregarded too. Question on the 
presence or absence of the quasiparticle branches in the 
K~ spectrum depends on the kaon energy and on the 
strength of the s- and p-wave attraction [14]. The short- 
range baryon-baryon correlations, which, as a rule, sup¬ 
press attraction, were disregarded in [6, 13] for simplifica¬ 
tion. A particular role of the correlations for the p-wave 
has been taken over in Ref. [15]. However, the relative 
strength of s- and p-wave attraction remained model de¬ 
pendent because no systematic investigation of the kaon- 
nucleon interaction including s- and p-waves has been 
available that time. 

Recently, the AT”-nucleon scattering has been stud¬ 
ied in the framework of the relativistic chiral SU(3) La- 
grangian imposing constraints from the A'“'"-nucleon and 
pion-nucleon sectors [16]. The covariant coupled-channel 
Bethe-Salpeter equation was solved with the interaction 
kernel truncated to the third chiral order including the 
terms which are leading in the large limit of QCD. 
All SU(3) symmetry-breaking effects are well under con¬ 
trol by combined chiral and large Nc expansions. This 
analysis gives an opportunity to extend results of [6, 15] 
taking into account off-pole (regular background) contri¬ 
butions to the kaon self-energy. The accurate fit of exper¬ 


imental data achieved in Ref. [16] fixes the values of the 
kaon-nucleon-hyperon coupling constants. Particularly, 
the S*(1385)-pole contribution to the kaon-nucleon scat¬ 
tering was proved to be sizable, being not included in 
works Ref. [6, 13]. 

The discussion of the s-wave and, especially, of the 
p-wave kaon-baryon interactions in nuclear matter is 
important for the kaon production in heavy ion colli¬ 
sions [17]. The momentum dependence of kaon yields 
is experimentally measured [18]. Also, the multi-branch 
K~ spectrum can be tested via ^-scattering on atomic 
nuclei [15]. Peculiarities of the AT^-nucleon interaction 
near mass shell are of great importance for the physics 
of K~ atoms [19, 20]. 

In this work we continue the study of the s- and p-wave 
Ar“-baryon interactions in dense baryonic matter of ar¬ 
bitrary isotopic composition. In particular, we present 
our results for the composition typical for the core of 
neutron stars. Our strategy to verify possibilities of s- 
and p-wave condensations is as follows. For given to¬ 
tal baryon density and the composition of the neutron 
star matter we find the K~ energy at lowest branch of 
the dispersion equation as a function of the momentum 
and then find the minimum in respect to momenta. At 
density, when w™'"(fcm) meets the electron chemical po¬ 
tential He, the reactions (1) become to be possible and 
the system undergos a second-order phase transition de¬ 
veloping a classical field with the K~ quantum numbers. 
If this state is realized for km = 0, we deal with the s- 
wave kaon condensation and, if km ^ 0, we deal with the 
p-wave condensation. 

Since solutions depend on many rather uncertain pa¬ 
rameters, we study possibilities of s- and p-wave conden¬ 
sations separately. For instance, discussing possibility of 
the s-wave condensation we put fc = 0 in the solution 
a'™“(fc), assuming that parameters of the p-wave inter¬ 
action are such that they do not allow for the p-wave 
condensation at the critical density for the s-wave con¬ 
densation. Then, trying to find a most appropriate pa¬ 
rameter choice we investigate how the possibility of the 
s-wave condensation depends on a parameter variation. 
Next, in the same manner we study a possibility of the 
p-wave condensation. Then we find the energy and pres¬ 
sure of the system with the kaon condensate at different 
spin-isospin compositions and densities. The system se¬ 
lects the spin-isospin composition corresponding to the 
minimum of the energy. We compare the pressure of the 
system with the kaon condensate with the pressure of 
the neutron star matter without condensate. If at some 
density the pressure and the nucleon chemical potential 
of the system with the kaon condensate coincide with 
those values of the neutron star matter (other density 
and isospin composition) without the condensate the sys¬ 
tem may come to the condensate state by the first order 
phase transition. Such a transition starts at a density 
described by the Maxwell construction. It occurs if an 
effective surface tension parameter is rather large. Oth¬ 
erwise, if an effective surface tension parameter is not 
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too large, a mixed phase is realized starting from even 
smaller densities. 

In this paper we describe the baryon matter in terms of 
the relativistic mean-field model (Sec. II). In Sec. Ill we 
introduce the kaon-nucleon interaction in vacuum, as it 
follows from the partial wave analysis of Ref. [16]. Then 
we separate the pole contributions of A(1116), S(1195) 
and S*(1385) hyperons in p-waves. Sections IV through 
VII are devoted to the construction of the kaon polar¬ 
ization operator. We start in Sec. IV with the polar¬ 
ization operator in the gas approximation but including 
the mean-field potentials acting on baryons. Besides the 
A(1116), 11(1195) and S*(1385)-nucleon-hole contribu¬ 
tions it contents a regular attractive part, weakly depen¬ 
dent on the kaon energy. In Sec. V we separate s- and 
p-wave parts of the kaon polarization operator. Occu¬ 
pation of hyperon Fermi seas is incorporated in Sec. VI. 
Repulsive baryon-baryon correlations are evaluated and 
included in the hyperon-nucleon particle-hole channels 
and in the regular part of the polarization operator in 
Sec. VII. In each section above we illustrate the strength 
of new terms included into the polarization operator and 
suggest effective parameterizations. We relegate the dis¬ 
cussion of contributions from kaon fluctuations (baryon 
self-energies, multi-loop corrections) to Appendix D. We 
argue that these effects do not modify substantially the 
kaon polarization operator at zero temperature in the re¬ 
gion of small kaon energies and momenta, which is of our 
interest here. In Sect. VIII we analyze different pos¬ 
sibilities of second- and Hrst-order phase transitions to 
the s- and p-wave K~ condensate states. In particular, 
we argue for the p-wave K~ condensation at p 2po 
(po — 0.17 fm“^ is density of nuclear saturation) arising 
via a hrst-order phase transition. In this phase tran¬ 
sition all the hyperon Fermi seas are melted and neu¬ 
tron star matter becomes proton-enriched (with approxi¬ 
mately symmetric isospin composition, N ~ Z). Depen¬ 
dence of the results on the specihcs of the EoS and the 
corresponding particle composition is illustrated in Ap¬ 
pendix A. Some technical information on the Green’s 
functions and the Lindhard’s function is deferred to Ap¬ 
pendices B and C. 

Throughout the paper we use units h = c = 1. 


II. BARYON INTERACTION IN 
RELATIVISTIC MEAN FIELD MODEL 

A. Lagrangian of the Model 

We consider a dense system consisting of baryons and 
leptons, which we describe by the Lagrangian density 
containing the baryon and the lepton contributions, re¬ 
spectively: C — Cb + Cl . 

The baryon matter at densities relevant for neutron 
star interiors is convenient to describe by the mean-held 


solution of the Lagrangian [21]: 

Cb = '^B{i ^ - g^B^ - gpB ^ -tB - ruB + gaB(y)B 
B 

df_,ad^a rnNKgaNcr)^ c{gaN<jf 

^ 2 2 3 4 

4 2 4 2 ’ ^ 

where all states of the baryon {J^ = i^) octet B = 
{n,p, A , , 5“°), interact via exchange of scalar, vec¬ 
tor and isovector mesons cr, p^. Heavier baryons 

do not appear at the baryon densities under consider¬ 
ation (p < 6po) and are not included, therefore. In 
(2) Ib denotes isospin operator acting on the baryon 
B. The held strength tensors for the vector mesons are 
given by = d^uju — ujp. for omega mesons and 
Piiu = dp, Pu — du Pfi for rho mesons. The equations of 
motion for baryons following from (2) give 

Eb{p) = eB{p)+ Vb , CBip) = +p‘^ , (3) 

where m*g = {iriB — g^B o') is the baryon effective mass, 
and Vb = gc,jBUJo + gpB Po 3 h b , with cr, ojq , pos being the 
mean held solutions of equations of motion for the meson 
helds. The composition of the cold neutron star matter 
at baryon density p is determined by the /3-equilibrium 
conditions Eb{pfb) = Pn — Qb Pe- Here pps is the Fermi 
momentum, gs is the electric charge of the given baryon 
species B, and p„ , pe are chemical potentials of neutrons 
and electrons hxed by the total baryon density p and the 
electroneutrality condition. 

The lepton Lagrangian density is the sum of the free 
electron and p~ meson contributions Ci = Ce + . 

The /3-equilibrium requires that the muon chemical po¬ 
tential equals to the electron one, p^ = pe, and, there¬ 
fore, muons appear in the system when pe exceeds the 
muon mass, . Electrons can be considered as ultra- 
relativistic particles. 

The energy density of the system is given by 


Et, 


Err 


E^^’^ipFB) 




Ernes + E^™(PFb) + Ei{pe), 

B /— 

^bniN (gcrNo)^ + 2<^{gaNo)^ 
lo 9 I99 I99 

+ 2 2 

PFB 

^ J dpp'^eBip), 

0 

_ 

■^9{pi-mi) J de ^- mf, (4) 

mi 


where 0{x) is the step function. 



4 


B. Coupling Constants 


The coupling constants are adjusted to reproduce 
properties of equilibrium nuclear matter: saturation den¬ 
sity po) binding energy -Ebind, compressibility modulus 
K, and the effective nucleon mass m’^{po). In the fol¬ 
lowing we take the value po = 0.17 fm“^, and the bind¬ 
ing energy iilbind = —16 MeV. For the nuclear compress¬ 
ibility modulus we take the value K =210 MeV, which 
follows from the variational calculation [22]. Following 
Ref. [21] we adopt the symmetry energy Osym = 36.8 MeV 
which lies well in the interval allowed by the microscopic 
calculations [23]. We take the effective nucleon mass 
= 0.85mAT, cf. argumentation in Ref. [24]. Of 
course, the values of the input parameters can be selected 
differently, in order to Ht experimental data. We refer at 
this point to the recent work [25] discussing the uncer¬ 
tainties. 

The corresponding coupling constants of the La- 
grangian (2) are: 


2 2 
9luN 

ml 
2 2 
gpN m-N 


= 54.6041. 


2 2 
9aNmM 

ml 


= 164.462, 


= 121.690, 5 = 0.0202832, 


c = 0.0471633. 


In order to verify the sensitivity of the results to the de¬ 
tails of the EoS we explore in Appendix A another set 
of the nuclear matter parameters, which we fit to re¬ 
produce the microscopic calculations of Urbana-Argonne 
group [26]. 

Including hyperons, one also has to specify hyperon 
couplings to the meson Helds. They are expressed via ra¬ 
tios to the nucleon couplings, xmh = 9mh/9mn with 
M = {a,u},p) and H = (A, E,5). Couplings to the 
vector mesons are estimated from the quark counting 
as = Xps = I and Xu^s = Xp's = |. Alterna¬ 

tively, relying on SU(3) symmetry [27] one would find 
XpT, = Xp’E = 1 for the p meson couplings. The scalar 
meson couplings can be constrained by making use of 
hyperon binding energies in infinite nuclear matter at 
saturation [28], extrapolated from hyper-nucleus data. 
Denoting the binding energy of a hyperon by we de¬ 
rive the following relation between the scalar and vector 
couplings: 

^Wnd = (glN Po/ml)Xu;H - {mN - m*j^)x„H 
= (80.73 - 140.70 MeV. 

There is convincing evidence from the systematic study 
of hypernuclei that for A particles — —30 MeV [29]. 
For S hyperons in nuclei, on the other hand, the data are 
still controversial, giving the broad band, —10 MeV < 
A^jnd < MeV, from a slight attraction to a strong 
repulsion [30]. Following Ref. [31] we adopt for 5 the 
values ~ —18 MeV advocated in Ref. [27]. 


To cover different possibilities we consider four cases: 

(I) vector meson couplings are taken according to the 

quark counting, a;(jA(E) = = 2 x^^s = 2 Xps = | and 

^^bind = -30 MeV, = -18 MeV, = 30 MeV; 

(II) the same as in case (I) but E^^^^ = —10 MeV; (III) 
the same as in case (I) but the rho meson couplings are 
taken according to SU(3) symmetry, i.e. Xps,s = 1; (IV) 
the same as in case (III), but for = —10 MeV. 


C. Particle Concentrations 

In Fig. 1 we show the resulting concentrations of dif¬ 
ferent baryon species, as function of the baryon density. 
Columns (I-IV) correspond to the four above choices of 
the hyperon-meson coupling constants. We see, that hy¬ 
perons appear in neutron star matter in all cases at den¬ 
sity p > pc,H — 2.5 -;- 3po- The latter critical value is 
rather insensitive to various choices of hyperon-nucleon 
interactions. However, the order, in which hyperons pop¬ 
ulate the Fermi seas, depends crucially on the details of 
hyperon-nucleon interactions. The hyperons do not 
appear, at least, up to 6po except for case II. However, 
even in this case their concentration is very small. The 
place of is readily taken by 5“ and A hyperons. This 
observation is in line with the results of work [32]. In 
case IV, we chose < 0, nevertheless the E hyperons 
do not appear due to increasing repulsion mediated by p 
meson with larger coupling constants than in case H. We 
see that in all cases I-IV the proton concentration and 
the electron chemical potential saturate with the filling of 
hyperon Fermi seas, in line with the results of previous 
works. We also see that all choices I - IV do not sup¬ 
port effects observed in [13], where A hyperons become 
more abundant than protons already at 3po- Therefore, 
the p-wave K^K~ condensation discussed in [13] does 
not show up in the framework of our model for all four 
parameter choices. 

Although concentration of hyperons is rather small in 
all the cases, their presence may have significant conse¬ 
quences affecting critical densities of p-wave K~ conden¬ 
sates. 


III. A -NUCLEON INTERACTION IN 
VACUUM 


The kaon-nucleon interaction in vacuum results as the 
solution of the coupled-channel Bethe-Salpeter equation 



where the triangle is the full meson-baryon scattering am¬ 
plitude with strangeness (-1) and the circle stands for the 
interaction kernel. This equation involves rescatterings 
through all possible meson-baryon intermediate states al¬ 
lowed by the strangeness conservation. The interaction 
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P/Po 

FIG. 1: Concentration of baryon species in neutron star matter. Panels (I-IV) correspond to four choices of the hyperon-nucleon 
interaction constants specified in text. 


kernel can be derived from SU(3) chiral Lagrangians [33- 
36] or phenomenologically adjusted to fit the data within 
the K-matrix formalism [37, 38]. 

In view of the general interest to the s-wave kaon con¬ 
densation prevailing in the literature so far, the most 
attention has been paid to the kaon-nucleon interaction 
in the s-wave. Although far below the threshold the s- 
wave kaon-nucleon scattering amplitude is rather smooth 
function of the kaon energy, close to the threshold the am¬ 
plitudes vary strongly due to the A(1405) resonance (see 
Fig. 2). Therefore, extrapolation of the scattering ampli¬ 
tude, adjusted to fit the data above the K~N threshold, 
into the subthreshold region crucially depends on the mi¬ 
croscopic model applied. Its in-medium modification is a 
matter of debate, too [39-43]. 

Up to recent, only a scarce information on the p-wave 
kaon-nucleon interaction was available. In the isospin- 
zero channel the small p-wave amplitudes were not sep¬ 
arated from the large contribution of the A(1405) res¬ 


onance in the s-wave, which dominates near threshold 
energies. In the isospin-one channel determination of the 
p-wave amplitudes remained also uncertain due to a lack 
of direct experimental information on the K~n scatter¬ 
ing at low energies. This gap has been filled by Ref. [16], 
where the relativistic chiral SU(3) Lagrangian imposed 
extra constrains from the iF^-nucleon and pion-nucleon 
sector. This analysis provides reliable estimates for both 
the s- and p-wave K~N scattering amplitudes, which we 
will use in the following. 


A. Forward Scattering Amplitudes 

The vacuum KN forward scattering amplitudes in the 
given isospin channel I have contributions from s- and 
p-partial waves: 


J 


=T®(s) -f T™(s) 

Q\s,p\k^) 


E{s,m%,mj^) +mN 
2niN 

{s — — Ak^ 

As 


F^\s) + Q\s,m%,m]^) {s) 


(I)/ 


( 6 ) 


r 


where s = {pA-k)"^, andp = {eN{p),p), k = (tu, fe ) are 4- 
momenta of the incoming nucleon and kaon, respectively. 
For nucleons and kaons being on mass-shell, the quan¬ 
tity is nothing else but the square of the 

center-of-mass momentum in the kaon-nucleon channel, 
E{s,m1j,mj^) = -|- Q'^{s, mj^, mj^) is the nucleon 

energy in the center-of-mass frame. 


rriN and rriK are the free nucleon and kaon masses. We 
shall neglect the isospin symmetry breaking within kaon 
and nucleon isospin multiplets, which are irrelevant in 
dense nuclear matter. 




s +p^ — k"^ 
2-\/s 


The invariant partial-wave amplitudes Fg'^ and Fp^ 
are related to the standard partial wave amplitudes 














6 


(cf. [44], section 3.1 for definitions) as 


Fi\s) 


8 7r^ 


mN 


■E{t 


‘AO '"‘K 


/o+(s)> 


8Try/s/Q‘^(s,m%,mj^) 
ruN + E{s,m%,m\) 


{41{s) 


■2/S(s)) . 


(7) 


The partial amplitudes fi± = in (7) are the 

scattering amplitudes for the given angular momentum 
I and total momentum J. Real parts of the partial am¬ 
plitudes t4 and t4 are shown in Fig. 2 and Fig. 3 by 
solid lines. Pronounced peak structures in the p-wave 
amplitude are due to the A(1116) pole in the isospin- 
zero channel and the E(1195) and E*(1385) poles in the 
isospin-one channel. 


B. Separation of the Pole Terms 


The hyperon s-channel exchanges are responsible for 
the most strong variation of the polarization operator at 
low frequencies and momenta. Therefore, they deserve a 
special consideration in view of the baryon modifications 
by the mean field-potentials. To treat them explicitly, we 
separate pole contribution from the p-wave amplitudes. 
Then we define the pole and pole-subtracted amplitudes 


5f4{s) 

— 


= -2- 
1 

Fllis) 

= -2- 
( 

4 2 {i 

3 s — 

m|. -1- ; 

The E* 

width. 

7^e(s) 4 

■ iKNis 

ttE and 

KN c 


[^(i) 

pole 

->2 
^KNA 




C'irMK 2 TOA 


(ruA -I- toat)^ 


{E{s,p'^, k^) + ttinY s ~ ’^A + * 0 ’ 

(ms + tunY 


( 8 ) 

(9) 


Fknt: 2 fns 


{E{a 


'-AT’ 


TTlN 


s — mi. + 


zO 


2C'|:^s* E{m^.,m%,ml^) + mN 


■7s*(s) E(s,m%,mj^)+mN 


.( 10 ) 


includes 


contributions 


of ttA . 


E(s,m‘^,m'^))\Q{s,mg,mQ)\^/{12TTy/s), where ms 
is the mass of the corresponding baryon B = 
(n,p, A, , S“°) and m^ is the mass of the corre¬ 
sponding light meson (j) = The coupling con¬ 

stants C',rA(E)E* can be extracted from the partial width 
of the S* hyperon. The values of the KNA{I]) coupling 
constants Ckna — —0.68/m^, Cknh — 0.34/m.n. enter¬ 
ing (9,10) follow directly from the amplitudes shown in 
Fig. 3. The values of Ckna and Ckny. are defined as 
couplings of A to N) and to isospin 

states, where AT, N stand for isospin-doublets, 11“ is 


the isospin-triplet, and r“ are the isospin Pauli matri¬ 
ces. Note that the value of Ckna, obtained on basis of 
Ref. [16], is smaller than that fitted by Jiilich group [45] 
and then used in Refs. [6, 12]. 

The scattering amplitude nearby the S* resonance can 
be only approximately described by the last term in (10), 
since the corresponding amplitude results from a multi¬ 
channel dynamics, which generates the energy dependent 
self-energy of the S* resonance to be included in the de¬ 
nominator. This self-energy, however, can be neglected 
at y/s relevant for calculations below. Ref. [16] gives the 
value for the E* resonance coupling Cknt.* — 0.84/m,r 
defined by the Lagrangian term CKNT.*^*fi^d^K^T°‘N . 

The pole subtracted amplitudes (s) are shown in 
Fig. 3 by dash lines. We see that the amplitudes in both 
isospin channels are smooth functions of in the sub¬ 
threshold region ./s < 1300 MeV of interest. 

Nevertheless, we would like to notice that procedure 
of the pole subtraction is not unambiguous. The forms 
(9,10), chosen here, will later exhibit themselves to be 
convenient, since they allow direct comparison with the 
analysis of Ref. [6, 13, 14]. 

Considering further the kaon-nucleon interaction in 
dense matter, we also need to take into account the mean- 
field potentials acting on baryons according to the La¬ 
grangian (2). 


IV. R- POLARIZATION OPERATOR FROM 
SCATTERING AMPLITUDES 

A. Gas Approximation and Baryonic Mean-Fields 

Our further aim is to construct the retarded K~ po¬ 
larization operator in baryonic matter n)^‘(a;, k) related 
to the propagator 

D-^^{u},k)=u}^-k^ 

The spectral function is given by AK(oJ,k) = 
—2Imi4A:(w, fc). The quasiparticle branches of the spec¬ 
trum appear in some energy-momentum region, if there 
the kaon width Tk = —2Imn)^*(a;, fc) is much smaller 
than any other typical energy scale. Then one can 
put Tk 0 in the kaon Green’s function and the 
quasiparticle branches are then determined by equation 
TleD~^{u, k) = 0. 

With the help of the scattering amplitudes in vac¬ 
uum we may extract the causal polarization operator 
np^(w,fc) related to the partial wave amplitudes F® 
and e4 in the gas approximation after integration over 
the nucleon occupation function: 
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FIG. 2: Real parts of the s-wave kaon-nucleon scattering amplitudes in isospin-zero and isospin-one channels [16]. 
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FIG. 3: Real parts of the p-wave kaon-nncleon scattering amplitudes Fp in isospin zero and isospin one channels (solid lines) [16]. 
The dash lines represent corresponding pole subtracted scattering amplitudes (8). 


Ig—wave(^i 


I s—wave (w, k) + Ip 

—wave (w, fe), 

Is—wave,p(^i k) Is—wave, 

f 2d^p E{s,mj^,k^) + rriN (1 


(2^)3 2v^ 


Ip—wave(^5 ^') — 


N^P 

wave,n(^: ^') 

2 d^p E{s, m'jf, fc^) -I- niN 


- np{p) + F^^\s) nn{p) I , 


Ip—wave ,^') d~ Ip— 


( 11 ) 


( 12 ) 


(2 7r)3 2y/m% + 


pz 


Q^(s,TO^,fc^) np{p) + F^^\s)nnip)^ild) 


where ni{p) are the nucleon Fermi occupations, i = 
{n,p), at zero temperature np(p) = 9(pY,p — |p|) and 
nn{p) = d{pF,n - IpI), S = (wF -bp2)2 _ (fc-|-p)2. 

There are simple relations between causal (“-,-”) and 


retarded (“R”) two-point functions, like Green functions 
and polarization operators, see eq. (B2) in Appendix B. 
For zero temperature causal and retarded two-point func¬ 
tions coincide for positive frequencies. For T ^ 0 their 
























real parts continue to coincide, whereas imaginary parts R and C, if it does not lead to ambiguity. 

are different, but can be interrelated. Bearing this in Exploiting decomposition (8) of the p-wave scattering 

mind, we will further suppress for brevity the subscripts amplitude we present 

_I 


Ip-wave(w, k) — fc) + fc), 

as the sum of the pole and regular parts. The pole part is generated by the hyperon exchanges 

k) = k) + k) + fc), 

TPoie. , .X ^2 f 2 d^P E{ml, m%, k^) - (wa + vin? ^ . 

J (2TrR 2 a/to^+p 2 s^-m\ + iO 

TPoie/ , \ r<2 f 2d^p -niN (mT.+niNf ^ ^ 

(--fc) = -Ckn^] + 

ypole / ^ f ^ ^ P , Tnj\^ j k ^ TTlS (^7 TV ? fc^) 

(u;, ) = J ^ ^ - ml. + ^ (s) 

The regular part, Ij^i'wave) can be expressed as 

C-wave(^, k) = i;-®wave(^, k) + k), 

including the part of integral (13) evaluated with dFp^ from (8), 


(14) 


{np{p) + 2n„(p)) 


(15) 

(16) 

(17) 

(18) 


C-wave(^,fc) = - 


2d^p E{s,m%,P) +mN 2 , 2 x 


(27r)3 2 a /mjf + p^ 


■ Q^{s,mj^,k^) 


X { i (s) + (s)) np (p) + (s) n„ (p) 


and the non-pole contributions from the hyperon exchanges 
'^Ip-wave(w> k) = (5Ia(w, k) + (5lE(a;, k ), 


i5Ia(w,A;) — —C](A!A J 
6l^{u!,k) = —Cj(A!A J 


2 d^p niA y/s — m?Ai + k^ (ruA + uin) 


(2 7r)3 2 Vi y^mjf + p^ Vs + uia 
2 d^p Toe Vs - m% + V (ms + totv) 
(2 7r)3 2 Vi v^ + WE 


■np{p), 


{np{p) -f 2n„(p)). 


(19) 

( 20 ) 
( 21 ) 

( 22 ) 


To obtain the last relation we used 

E[s,niAf,k ) — E[nipi jITIa! V ) 


rriHVs - mV + V 
2mpfVs 


{Vs-mn), H = {A,'E}. 


(23) 


The above construction of the polarization operator, cor¬ 
responding to the gas approximation, does not take into 
account mean-field potentials acting on baryons, vertex 
corrections due to the baryon-baryon correlations, and 
possible modification of the scattering amplitudes in the 
medium. The modification of the baryon propagator on 
the mean-field level can be easily incorporated in the in¬ 
tegrals (12,13) by the replacement ttin m^- Effects 
induced by this modification in the kinematic prefactors 
in (12,13) can be easily traced, as we will demonstrate be¬ 


low. Scaling of nucleon mass in s is more subtle. Solving 
the coupled-channel Bethe-Salpeter equation one sums 
all the two-particle reducible diagrams for the part of 
the s-plane corresponding to K~N scattering. This ap¬ 
proach is explicitly crossing non-invariant and continu¬ 
ation of amplitudes far below K~N threshold can gen¬ 
erate artificial singularities in the scattering amplitude. 
In Ref. [16], from where we borrow the amplitudes, the 
approximation scheme for solution of the Bethe-Salpter 
equation, was furnished in such a way that the K~N 
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and K'^N scattering amplitudes exhibit the approximate 
crossing symmetry, smoothly matching each other for 
y/s ~ mN- Therefore, amplitudes depicted in Fig. 2 and 
Fig. 3 are still physically well constrained in the corre¬ 
sponding intervals of ^/s shown there. However for some¬ 
what smaller ^/s the K~N s-wave scattering amplitude 
gets unphysical poles. To cure this problem the complete 
solution of the Bethe-Salpeter equation for K~N scatter¬ 
ing has to be redone with the medium modified baryon 
masses. Fortunately, there are some indications that it 
would not change the results drastically. For the integral 
with the s-wave scattering amplitude we will demonstrate 
that the final results can be nicely modeled with the po¬ 
larization operator following from the leading-order chi¬ 
ral Lagrangian, which has now explicit dependence on 
the baryon masses. The loop corrections due to iteration 
of the interaction kernel should be suppressed for small 
kaon frequencies to keep approximate crossing invariance 
of the amplitude. The pole subtracted p-wave ampli¬ 
tude is rather smooth function of ^/s, as it is shown in 
Fig. 3, being mainly determined by contact terms of the 
chiral Lagrangian and, thus, has a weak baryon effective 
mass dependence. Hence, extrapolating the amplitude to 
somewhat smaller ^/s, we do not expect its strong vari¬ 
ation. Opposite, the part of the polarization operator 
generated by the hyperon poles, depends strongly 

on baryonic mean fields changing the pole positions in the 
amplitude. Therefore, this part will be treated explicitly 
in the course of our consideration. 


B. Pole Part of Polarization Operator 

Here, we find contributions to the K~ polarization op¬ 
erator from the hyperon poles in the K~ N scattering am¬ 


plitude determined by (15,16,17). Relying on explicit cal¬ 
culations of Refs [6, 15] we can easily incorporate scalar 
and vector mean-fields acting on baryons. The scalar 
field is taken into account with the help of the replace¬ 
ment rriB rn*g and the baryon vector potentials are 
included in the pole terms by the shift of the kaon fre¬ 
quency OJ ^ OJ + 6viH, with SviH = Vi- Vh, i = {n,p}, 
H = {A, E} , that immediately follows from the differ¬ 
ence of the baryon energies, see (3). Please notice that 
this energy shift is obvious only for the pole contribu¬ 
tion to the polarization operator. Generally, due to the 
absence of the gauge invariance for the massive vector 
fields such a shift is not motivated for more complicated 
diagrams. 

Writing down explicitly all contributions we cast 

n(pole.0)(^^^) = ipol^ave(u;,fc) 

= (c., fc) + njpf(c., fc) + 2 (c., k) 

(a;, fe) + 2 (c., fc), (24) 

where each term is equivalent to the pole contribution 
of the hyperon particle - nucleon hole loop diagram 
(Schrodinger picture) 


written in terms of the Lindhard’s function [74] as 



niff ’\^y^) — ^KNH [(m^ — m^)^ — (w-|-(5viff)^-f fe^] <I>iff (w, fc), 

2 / + m%)'^ - {u + Svin)'^ + k^ 

(w, k) = —i-• 


(26) 


We reserved the superscript ”(0)” for each term in (24) 
indicating that neither further self-energy of baryons be¬ 
yond the mean field nor the vertex corrections due to the 
baryon-baryon correlations are included. The (retarded) 
Lindhard’s function <I> accounting the relativistic kine¬ 
matics is defined as 


centrate further 


$iff(a;,fc) 

A)^ff(a;,fc,p) 


^ PPi ^ , 

f dpp A+^{uj,k,p) 

2Tr^\k\ J ei{p) ^ k,p)\ ’ 

[oj -I- 6viH + ei(p)] ^ — Cff (IpI T |fc|) )(28) 


$iff (w,fc) 


2d^p 


4m^ 


(2 7r)3 2 ei{p) s — rrijj -I- f 0 


riiip). (27) 


For the case of zero temperature, on which we will con- 


where ppi is the Fermi momentum for the given nu¬ 
cleon species. Non-relativistic form of the Lindhard’s 
function used e.g. in Ref. [24] (in different normaliza¬ 
tion), is recovered with the help of expansion eff(p) « 
m)j -|-p^/(2m)j) in (28). 
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The imaginary part of the (retarded) Lindhard’s func- Here are the upper and lower borders of the corre- 

tion is obtained as an analytical continuation ln(a;) = sponding particle-hole continuum 

In \x\ iiT 0{—x) leading to a non-zero contribution for 

< w < (29) 


J 




f \J{rn*fj - + 5viH 

y EH{\k\ +PFi) — Ei{pFi) 
EH{\k\ — PFi) — Ei{pFi). 




Baryon energies include vector potentials according to (3). Properties of Im<I> are discussed in Appendix C. 
An approximate expression for renders 


(30) 

(31) 


$i_f/(w,fc) « -- 






m*j^ 





fc 3 eFt 

ill 

2 


4 

5 

'uj,k,PFi), Ef 

i = ei{PFi), 



(32) 


being valid for to < The accuracy of this 

approximation is illustrated by Fig. 4. 

For completeness we also quote here the low momen¬ 
tum limit of (32), 


2 Pi m*^ 

■ k^ 

1 + — 

fl-b - 

)1 

epi ^iH 


V 5 



Re$i_f/(w,fc) 


= A)*)j(u;, 0,pf*), |fe| «C |Aj//(u;,0,pFj)|/(2pFi)(33) 


namely, 

ns^(w) = Is-wave(w, 0) + Ip_wave(w, 0), (34) 

Hp (cu. A;) = [Is_.,vave (^7 ^) Is—wave (^5 0 )] 

+ [Ip —wave(‘[^5^) Ip—wave (a;,0)]. (35) 

Next two subsections are devoted to discussion of the 
s-and p-wave parts of the polarization operator. 


where pi is the density of the given nucleon species n or 
P- 


A. P-wave Part 


V. S- AND P-WAVE PARTS OF THE 
POLARIZATION OPERATOR 


In our discussion we would like to put a particular 
emphasis to in-medium effects, which modify the K~ 
spectrum at finite momenta. For this purpose we dehne 
the momentum independent part, called s-wave part of 
the polarization operator, and the momentum dependent 
part, called p-wave part of the polarization operator, 


n^°)(c.,fc) 


n^°^(a;)-bn^°^(a;,fe) 


n[?)(o;,0) 




The term n[?^(a;,0) does not depend on fc, whereas the 

depends on fc, vanishing at 


term j^H^^ (w, fc) — np''(w, 0) 

|fc| = 0. In order to avoid misunderstanding we point out 
that the s- and p-wave scattering amplitudes contribute 
to both the parts (12), (13) of the polarization operator, 


( 0 )/ 


Following decomposition (14) we split the p-wave part 
of the polarization operator into the pole and the regular 
contributions 

n® (w, k) = k) + k). (36) 

For the pole p-wave part we have 

(cu, fc) = (cu, fc) - (w, 0) 

= (c., fc) + (cu, fc) + 2 (cu, k) 

+n(P’°)(a;,fc) + 2nf^?l(u;,fc), (37) 

(u;, fc) = (cu, fc) - (cu, O), (38) 

where we used (15)-(17) and (24). The contribution 
of S*-nucleon-hole term is rather sensitive to the val¬ 
ues of the mean-field potentials acting on E* which are 
unknown. Therefore, we investigate two choices in fur¬ 
ther. First, we assume that E* couples to the mean field 
with the same strength as the E hyperon (Vs. = Vs, 
5(tE* = (/cte)- In the second case, we detached E* from 
the mean field potentials {Vs- = 0, g^s- = 0). 
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FIG. 4: The real part of the relativistic Lindhard’s function (28) is shown by solid lines as a function of energy for different 
momenta. Dashed lines correspond to the approximate relation (32). The curves labeled with (1) are calculated for the Fermi 
momentum pp = 300 MeV and those labeled with (2) stand for pp = 400 MeV. 


Expanding the p-wave pole part of the polarization operator in small kaon momenta we have 

I 


= Vnh 0) + vlfH {(''tT'H ~ — (^ + ^VinY) -^ 

‘hjs* fe) = C’Itate* 0) 0)- 


(39) 

(40) 


The regular part of the p-wave polarization operator 
is defined by 


iip 


(w,fe) 


— [Is—wave(^i Is—wave(^; 0)] 

+ [lp‘l®wave(^,fe)-lp-®wave(^,0)]. 


At small kaon momenta the real part of can be 

written as 

ReU^P^^’°\uj,k) = (bp{uj)^ + bn{oj)^)+0{k^), 
\ On Pq/ 


PO 

h{uj) = bf\uj) + bf\uj), 


bf\uj) = — ^::^Is_wave,z(w,fc) 


Po d 


bT\u:) = 


P^ dk^ 
Po 

Pi dk 


2 ^p—wave,i(^i 


fc=0 


fc=0 


(41) 

(42) 

(43) 


where we used that the real part of the kaon polarization 
operator is even function of the kaon momentum. The 
quantities and are shown in Fig. 5 (left plane) 
for T = 0 as functions of the kaon energy for several val¬ 
ues of the densities pi. In these calculations integrals Is 


and Ip have been evaluated with the free nucleon masses. 
We see that these coefficients are almost density indepen¬ 
dent and only weakly dependent on the kaon energy in 
the interval 100 MeV < uj < 250 MeV of our interest. As 
we have discussed in the beginning of this section, we re¬ 
place the baryon masses by the effective masses, as they 
follow from the mean field solutions, only in the kinemat- 
ical prefactors in (12,13). The results are shown in the 
right plane of Fig. 5. The coefficients bn ’ ’ depend mod¬ 
erately on the density, as before, whereas the coefficients 
bp ’ exhibit a stronger density dependence, which can 
be parameterized by the factor rriN/m’^, as it is demon¬ 
strated in Fig. 5. Energy dependence remains to be weak 
in the interval 100 MeV < uj < 250 MeV of our interest. 

Our result (41) is derived for rather small values of 
kaon momenta, | fc |<C rriK- In order to find the ac¬ 
tual value of the p-wave K~ condensate amplitude in 
most general case one needs to deal with momenta up to 
I k 1^ pY^n ^ rriK- To satisfy the latter general case we 
extrapolate our result for the regular part of the polar¬ 
ization operator to such momenta. Luckily, within our 
approach critical points of s- and p-wave condensations 
deviate not as much from each other, as we will show it, 
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FIG. 5: Left panel: Coefficients (42) and (43) of expansion (41) of the regular p-wave part of the polarization operator. Solid, 
dash and dotted lines show calculations done for nucleon densities Pn,p = 1,3, and 5po, respectively. Right panel: the same as 

(S P) 

in the left plane but integrals are evaluated with the effective nucleon mass in kinematical prefactors and coefficients bp ’ ' are 
scaled by m)(,/mjv. 


and in this case the kaon condensate momentum in the 
vicinity of the critical density has extra smallness. Also, 
the main contribution to the kaon polarization operator 
comes from the pole terms, which are written explicitly 
for arbitrary momenta (24). Thereby, ambiguity of men¬ 
tioned interpolation should not significantly affect our 
conclusions. 


B. S-wave Part 

The kaon-nucleon interaction determines the following 
contributions to the s-wave part of the K~ meson polar¬ 


ization operator, 


Hg ^(w) — ls-wave(w,0) -h 0) 

(w, 0), (44) 


where the last two terms correspond to non-pole and pole 
parts of the hyperon exchange terms in the amplitude, 
respectively. The term is given by (24). Using 

(20) we present as follows 






0) 




(w,0) -h2m 




with 


^(w,o) — 


2d®p {m*jj + m*j^Y , x 

(2 7r)3 2ei(p)A/s(l(v^ + ’^ff) ^ 


(45) 


where sp = (w -|- ei(p))^ — and we also included dependence of effective masses on the mean field. For T = 0 and 
for small kaon energies, the integral (45) can be well approximated by the following expression 




(reg,0) 

iH 


(^,0) 


'-'KNH 


( - rn-N 

V 


seal 

Pi 


+ UJ Pi 


CO 




‘AT 



(46) 


where stands for the scalar density of nucleons defined by 


f 2d®p m% 

J (2 7r)3ej(p)' 


0 


(47) 
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According to (26), (28), the last term in (44) can be cast as 


n(p°Ho)(^,o) = nJ,p°^'=’°^(w,o) + nJ,P2°^'=’°^(w,o) + 
niff ’°^(‘^)0) = C'l^ATff ((wff ~ ■'^ff)^ ~ (^ + 

(m^ + m^)" r dpp"^ 

27r2 J eM[^U^,0,P)+^0] 

that follows from (28) at |fc| ^ 0. The imaginary part, lmnl^°^*^’°^(a;,0), is given by 


ImnlP°‘®’°^(w, 0 ) = -iC^NH - {uj + SviH)^) 

(m^ + mj^)2 


87r(a; + Sv^h)'^ 


‘ff 


-(r 




6u. 


iH) 


(48) 


(49) 


being non-zero for energies Wjj:^(0) = Eh( pFi)—EiiPFi) < 
uj < co+fj{Q) = Eh{0) - E,{Q). 


C. Energy of the Lowest Branch of the Dispersion 
Equation at fc = 0 

In this sub-section we illustrate strength of different 
terms in (44) applying it to the problem of the s-wave 
kaon condensation. 

The neutron star matter becomes unstable with re¬ 
spect to reactions (1) with the production of the zero- 
momentum K~ meson, when solution (ws = w™“(fc) at 
fc = 0) of the dispersion equation 

“ Rens(ws) = 0 (50) 


cf. [8], 


^S^'^iPn,Pp) = \jmH - Sk + V^ - Vk , 

c ^ / V' / seal I scal\ , /-i i seal 

[EkN [Pp + Pn )+C (Pp - 
Tr _ (2 Pp + Pn) 

. _f 9 , 



(54) 


where / ~ 90 MeV is the pion decay constant in the chi¬ 
ral limit [16], and Ekn and C stand for the isoscalar and 
isovector kaon-nucleon S-terms related to the explicit 
chiral symmetry breaking. The SU(3) symmetry predicts 
C = m\ (2toh — 3mi;-|-mA)/(16 [rn^H — rn^)) « 66 MeV, 
cf. [4]. The model polarization operator leading to the 
dispersion relation (54) can be written as, cf. [8], 


related to the lowest branch of the spectrum, meets the 
electron chemical potential. Then the s-wave K~ conden¬ 
sation may occur by the second-order phase transition. 

In Fig. 6 we present the energy at the lowest K~ 
branch of the dispersion equation (50), as function of 
the density at the momentum fc = 0 [75]. The hyperon 
interaction is taken according to case I. In order to illus¬ 
trate the strength of different contributions to the s-wave 
part of the polarization operator (44) we consider several 
test polarization operators 

n('’°)(a;) = I„e(ui,0)+I^'L®wave(t^,0), (51) 

^ n^^’°)(w)+OT(‘'"S’°)(w), (52) 

n®(a;) = n^^’°^(w)-hn(p°'®’°)(w,o). (53) 

In panel (a) of Fig. 6 solid lines show the energy of the 
lowest branch solution of (50) with Ilg^’^^ for the cases 
of pure proton and neutron matter. The contribution of 
Ip®®wave(^j 0) is found to be very small, at the level of few 
percent. It is instructive to compare this our result with 
that given by the frequently used parameterization of the 
K~ spectrum motivated by the leading-order chiral per¬ 
turbation theory (yPT) expansion of K~ N interaction. 


n: 


(XPT.O) 


(w) = - 


^KN ^ =eal 


^scal 


c 


P 

Pp + p 


■ LO . 


) - p ^Pl - pT) 

(55) 


Spectrum (54), calculated using Ekn = 150 MeV, is 
shown in Fig. 6 (panel (a)) by the dash lines. We ob¬ 
serve a good agreement of the model spectrum with 
that follows from numerical evaluation of the integrals 


F 


,(u;,0)+l^is_,,(u;,0) 


. Note, that obtained value 


of the effective kaon-nucleon S-term is 2-3 times smaller 
than that used as an ad-hoc parameter in Ref. [8], where 
the same parameterization (55) has been exploited. 

The results for the realistic composition of neutron star 
matter, shown in Fig. 1, are presented in panel (b) of 
Fig. 6 for case I of the hyperon-nucleon interaction as a 
representative example. The solid lines depict the energy 
at the lowest branch of the dispersion equation calculated 
for fc = 0 with Ilg^’®^ and Hg^^ Dash lines show 

solutions obtained with the approximate expressions (55) 
in Ilg^’®^ and (46) for (5nF®s>o) Hg^’^^. Excellent coin¬ 

cidence of the curves justifies the accuracy of (55) and 
(46). 
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p/p„ p/p„ 


FIG. 6: Panel (a): The energy of the lowest K~ branch of the dispersion equation (50) at fe = 0 calculated with Ilg^’^^ (solid 
lines) for proton and neutron matter. Dash lines present approximate spectra given by (54). Plane (b): The energy of the 
lowest K~ branch of (50) at fe = 0 for real neutron star matter (case I of the hyperon-nucleon interaction) for the s-wave 
polarization operators Ilg^^ given by (51,52,53) (solid lines). Dash lines are solutions obtained with approximate 

relations (46,55). Dotted line shows the electron chemical potential. Dash-dotted line depicts the border of the imaginary part 
of the K~ polarization operator (border of hyperonization). Short-dashes show the boarders of the regions of non-physical 
solutions of (50) characterizing by negative sign residues. 


Crossing point of the solid and dotted lines corresponds 
to the critical density of the s-wave condensation for the 
case of a realistic neutron star composition. We observe 
that the lines corresponding to do not meet the 

chemical potential (dotted line). Therefore, the s-wave 
kaon-nucleon interaction, following from Ref. [16], would 
not support a second-order phase transition into the s- 
wave K~ condensate state due to the small value of the 
kaon-nucleon sigma term following from the analysis [16]. 
However, an additional attraction comes from the term 
'^Ip-wave included in It makes the reactions (1) 

possible at density > 4.5po- Another attractive piece is 
the pole term Ip[]]^j^.^,g(a;, 0) taken into account in H®. 
The significance of these terms, originating both from 
the hyperon exchange diagram in KN interaction, was 
pointed in Ref. [6]. Both mentioned contributions were, 
however, disregarded in works [5, 7-12] discussing s-wave 
K~ condensation. 

The curves cus calculated with the full s-wave polariza¬ 
tion operator H® have cuts. In the region between the 
cuts equation (50) has no solutions with positive residues, 
cf. [6]. The dash-dotted line depicts the border of the 
imaginary part of the K~ polarization operator. We see 
that, fortunately, within our approach the curve calcu¬ 
lated with Hg®^ and /ig meet at energy below the region 
of the imaginary part. Thus, with the full polarization 
operator (44), (53) we recovered statement of previous 
works [5, 7-12] (where, however, the 2 4-3 times larger E 
term was used) on the possibility of the K~ condensate 
production in reaction (1) at rather moderate densities, 
Pc S — 2.7po in our case. The reader should bear in mind 


that baryon-baryon correlations are still not included in 
the above analysis. They will increase pc.s- This issue 
will be addressed in section VII. 


VI. CONTRIBUTIONS OF THE HYPERON 
FERMI-SEAS TO THE POLARIZATION 
OPERATOR 

When nucleon density exceeds the critical density of 
hyperonization Pc,H^ tbe Fermi sea of the hyperon H be¬ 
gins to grow, and the K~ polarization operator receives 
new contributions. 


A. Pole Terms 

New contributions to the K~ polarization operator re¬ 
late to the diagrams (in Schrodinger picture) 



where the hyperon and the nucleon interchange their 
roles compared to the above discussed the hyperon- 
nucleon-hole NH terms. At finite temperatures these 
diagrams contribute for all densities. But this contribu- 
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tion is exponentially suppressed for p < pc^H and small 
temperatures. 

Account of the hyperon contribution to the pole part 
of the polarization operator is simply done with 

the help of the replacement in (26) 

k) k) + w, —k) (57) 

where the last term implies interchange of all indices 
i ^ H in (27), (28), (32). This replacement results in 
appearance of an extra term 

+ 2n^^°y\-uj,-k), (58) 

to be added to the total polarization operator. This term 
depends on the hyperon density and contributes only at 
densities p > Pc,H- 


B. Regular Terms 



FIG. 7: The energy of the lowest branch of the dispersion 
equation at fe = 0 calculated with the polarization operator 
ns°hyp given by (61) is depicted by solid lines. Dashed line 
shows the corresponding solution with the polarization op¬ 
erator (53). The dotted line depicts the electron chemical 
potential. 


There are no any experimental constraints on the hy¬ 
peron contribution to the regular part of the polarization 
operator so far. As a rough estimation we may suggest 
an extension of the model polarization operator (55) to 
the hyperon sector according to the leading order terms 
of the chiral Lagrangian 

C 1 

/ ■*■ seal seal . seal\ 

- -pi-^PA - Py.- +Ps- ) 

I ^A seal I P'S- +2pH- , /rqx 

+ J2 Pa H- pp - 

In this expression we utilize the value of "Skn from the fit 
with the formula (54) to the numerical results in Fig. 6 
{T,kn — 150 MeV), whereas the values of coefficients 
C « 66 MeV and Ca = m|. (toh— mA)/(12 (m|.—m^)) « 
34 MeV and / ~ 90 MeV are predicted by the chiral 
SU(3) symmetry. To estimate the non-pole contribu¬ 
tion from the nucleon u-channel exchange (analogous to 
) we use approximate relations (46) 


(reg.O)/, A_ /^2 

yp 1 

*2 *2 


= -^kna 


m‘^ — m 


2m\ 




A ^scal 

Pa 


2 TO? 


-2C 


''S ^scal 

— Ps- 


■LO PK 

KNT. 

^ Pt.- 


TTls 


m 


■PA 


N 


ml 


■ m 


• Pt.- 


N 


(60) 


C. Energy of the Lowest Branch of the Dispersion 
Equation at fe = 0 

Solutions related to the lowest branch of the dispersion 
equation (50) for fc = 0 calculated with the polarization 
operator 

<hyp(^) = + 

+ m(P;;'=’°)(a;,0) (61) 

are shown in Fig. 7 by solid lines in comparison with cor¬ 
responding solutions obtained without inclusion of hyper- 
ons (dash lines). Calculations are done for hyperon cou¬ 
pling constants corresponding to case I. We see that pres¬ 
ence of hyperons produces an additional small attraction 
only at rather high densities (> 4po)- The reasons are 
partial cancellation of the attractive term and 

the repulsive term and that in the framework 

of our model for description of the neutron star matter 
hyperon concentrations are much smaller than the neu¬ 
tron concentration and even smaller than the proton one. 
Thus, population of the hyperon Fermi seas only slightly 
affects the s-wave part of the polarization operator. 

This allows us not to care much about the hyperon 
Fermi sea occupations considering fc = 0 case. 


VII. BARYON-BARYON CORRELATIONS 


With this estimation we do not take into account con¬ 
tributions to the regular p-wave part of the polarization Operating with the polarization operator constructed 
operator cx p//. by integrating the meson-nucleon scattering amplitude 













16 


over the nucleon Fermi sea, e.g., as in (12,13), one as¬ 
sumes that all multiple meson-nucleon interactions are 
independent from each other and have the same probabil¬ 
ity proportional to the nucleon local density p{r ). How¬ 
ever the successive meson-nucleon scatterings in dense 
nuclear matter are not independent because of the core 
of nucleon-nucleon interactions and the Pauli exclusion 
principle [47, 48]. The probability to find two nucleons i 
and i' at the positions ri and r 2 , respectively, is propor¬ 
tional to the two-particle density 

Pii'{ri.,r2) = [1 + - r2 \)] Pi{ri) Pi'{r2 ) 

with the correlation function Cu' (r) < 0 and is, there¬ 
fore, reduced in comparison to the product of two single 
particle densities. The correlation function can be ap¬ 
proximately written as 

(r) « + 5u' Cf [1 + (62) 

with contributions from the hard core, C'™’’®, and the 
Pauli exclusion principle, assuming that both cor¬ 

relations contribute multiplicatively. The former can be 
taken from the description of nuclear matter with the re¬ 
alistic nucleon-nucleon interaction. The convenient pa¬ 
rameterization was suggested in Ref. [49], C'®°''®(r) « 
—joirnor) with toq « where ji{x) stands for 

the spherical Bessel function. For the Pauli correla¬ 
tion we use expression for the ideal fermion gas [50], 
(^Pauii(^) ^ -9jj{pF^r)/{2pl.r^). 


In diagrams, correlation processes can be presented by 
symbolic equation 


B 


B 


B 


B' 




(63) 


K- 


where the wavy line relates to the kaon, the sum goes 
over the baryon species. The line without arrow means 
that both particles and holes are treated on equal footing 
(the conservation of charges, e.g. strangeness, baryonic 
number etc., in each vertex is implied). The hatched tri¬ 
angle is the bare scattering amplitude (scattering on the 
particle or on the hole) and the full triangle stands for 
the amplitude including baryon-baryon correlations. The 
dotted line symbolically depicts the two-baryon correla¬ 
tion function Cbb' due to the BB' correlations through 
the core and the Pauli principle. There are nor exper¬ 
imental information neither theoretical estimations for 
the hyperon-nucleon and hyperon-hyperon correlations. 
Since the latter ones are less relevant for our discussion 
below we will neglect them. Thus we actually include 
only minimal correlations given by Cu/ nucleon-nucleon 
correlation functions. 

We consider first the correlation corrections to the s- 
wave part of the kaon polarization operator, Hg^^ given 
by (44), and the regular p-wave part , from (41). 

We separate contributions induced by the scattering on 
a nucleon of a given isospin species i = {n,p}. 


A. Correction of s-wave and Regular p-wave Terms 

General derivation of the corrections to the meson 
propagation in dense nuclear matter due to the nucleon- 
nucleon correlations (so-called Ericson-Ericson-Lorentz- 
Lorenz corrections) can be found in Ref. [51, 52] for pi- 
ons. In Ref. [53, 54] it was extended for kaons. 


n(°)(o;) 


Is—wave, i (w, 0) -f Ip —wave,i (t^,0). 


r(0), 


cf. (12,13). Then, adopting results of Refs. [54, 55] we 
may present the polarization operator terms corrected by 
baryon-baryon correlations as 


ns(w) = 


n; 3 ®®(a;,fc) =fc' 


ns,n(w) + ns,p(a;) + 2nSj„(a;) nSjp(a;) (oj) 

l-fisMnsJcu) (4n\i^))^ 

2 bp (uj) + bn(u}) + 2bp(uj)bn(co)^ip\u;) 

1 -bp(u;)bn(uj)(^ip\uj))^ 


ns.i(w) = 


bi(aj) = 


nS(^) 


(Pi/po) 


l-bi{uj)p^^lP{uj)/po 


(64) 

(65) 


Functions, and are defined as, cf. Ref. [54], 

= J d^rCii,{r)D°j^{uj,r), ( 66 ) 

= IJ d^rCiAr)V^ D%{u;,r) 

= -^ (Cii'(O) - (ml: - ^lv^(w)) , (67) 


containing Pauli and core contributions fol¬ 
lowing (62), i,i' = {n,p}, and D^{uj,r) = 
— exp(—i/ml- — r)/(47rr), as the free kaon propaga¬ 

tor in the ’’mixed” representation. 


Using (63) one finds that the repulsive core contributes 



with 


(P, core) _ ^(P,core)^^ _ 1 




K 


^) = o 


3 {ml + mj^ — 


( 68 ) 


to the p-wave ^’s and with 


If we approximate the free KN scattering amplitude 
(hatched triangle) in (63) by the hyperon-exchange di¬ 
agram, the same one, which produces the particle-hole 
diagrams (25), we can see that the account of correla¬ 
tions via (63) is equivalent to the replacement [76] 


^(S,core) _ j^(S,core) 


(w) = 


{ml + 


(69) 


to the s-wave ones. The contribution from Pauli spin- 
correlation (second term in Cw ) to the p-wave correlation 
function is shown in Fig. 8 (left panel) as a function 
of the Fermi momentum for different kaon energies. We 
see that the correlation parameter decreases with density 
since the baryon-baryon correlations hold baryons apart 
from each other and suppress, thereby, effect of the Pauli 
exclusion principle. The right panel of Fig. 8 presents 
the values of the correlation function (66, 67) calculated 
for w = /ie in the neutron star matter with the hyperon 
coupling constants according to case I. 

We leave the contributions from the hyperon Fermi 
seas to the regular part of the polarization operator, 
(5ng°|^yp, without corrections due to the baryon-baryon 
correlations in view of their small contributions. 


B. Correction of p-wave Pole Terms 

We turn now to consideration of effect of correlations 
in the particle-hole channel. 



with a modified vertex (fat point) obeying equation 



where the particle-hole irreducible box Tj^^, (light 
square) can be expressed in terms of and kaon- 
nucleon-hyperon coupling constants. 

Below we would like to put the discussion on a more 
phenomenological level. According to the argumenta¬ 
tion of the Fermi-liquid theory [56], the particle-hole ir¬ 
reducible box ^ weak dependence on incoming 

energies and momenta and can be, therefore, parame¬ 
terized in terms of phenomenological Landau-Migdal pa¬ 
rameters: 
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= Coff^, 


(72) 


with a = 1,2,3. The amplitudes are normalized with 
Co = 300 MeV • fm^, that allows to compare the val¬ 
ues for the hyperon-nucleon correlation parameters with 
those for the nucleon-nucleon correlations introduced in 
Ref. [24]. In (72) P^^^^ = (1 — fs • t)/ 3 is the projec¬ 
tion operator onto the given UN state with isospin 1/2, 
ts are the isospin-1 matrices and r are the Pauli ma¬ 
trices of the nucleon isospin, as above. The projector 


Pe*p^^ is defined analogously. The spin-spin operators in 
the particle-hole channel are given by Si 2 = (cricr 2 ), and 
<S'i 2 = (S^i>S' 2 ^), with S' standing for spin-operator, which 
couples spin 1/2 and 3/2 states. 

The inclusion of correlations according to (71) brings 
the pole polarization operator (24) into the form 
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FIG. 8: Left panel: contribution of the Pauli spin-correlations to the p-wave correlation function (67) vs. the Fermi momentum 
for different kaon energies, is the Fermi momentum for p — po- Right Panel: the s-wave (dash lines) and p-wave (solid 
lines) correlation functions (66,67) evaluated for ui = pe in neutron star matter (case I) as function of the total baryon density. 
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(73) 


C. Correlation Parameters 


This can be shown in diagrams as 


To our best knowledge there is no direct experimental 
information on the values of the Landau-Migdal parame¬ 
ters for the hyperon-nucleon interactions /)(, 52 (e»)i and 
/as- In principle, this information could be extracted 
from the data on multi-strange hyper-nuclei, which, how¬ 
ever, are rather poor, if not absent. In the work [15] 
the Landau-Migdal parameter /)( was estimated in line 
with Refs. [57], where Landau-Migdal parameters of 
the nucleon-nucleon interaction were calculated within 
Ericson-Ericson-Lorentz-Lorenz approach. We will fol¬ 
low this approach, estimating these parameters. Further 
corrections can be computed, as in Ref. [49]. 

Following [57] we assume that the squared block in (72) 
is determined by exchanges of the kaon and the heavy 
strange vector meson K* with the mass rriK- — 892 MeV. 



Intermediate states in the processes (74), (63) involve 
large momenta suppressing medium effects. In this ap¬ 
proximation, correlation parameters are equal for HN~^ 
and NH~^, for the nucleon and the hyperon of given 
species. Then, including hyperon Fermi seas, the pole 
term of the polarization operator is corrected with the 
help of the replacement (57) in (73). The block (74), 
being evaluated at zero momentum and energy transfer. 
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contributes to the local interaction in (72) as 


rp\oC ^(PTCOre)/-N 

Ihn ~ ^KNH ^KNH 1^1 

+ Ck*NH Ck-NH ^ 

i7 = A,E,E*, 

^a°e(^) ~ CknaCkn'^Ck’^°’^‘^\^) 

+ Ck*naCk*n’-s^k*'^°^^\^) ■ C ^ 6 ) 


For shortness we do not write here explicitly the spin 
and isospin operators which are exactly the same as in 
(72). The vector meson coupling constants Ck-nh in 
(75) correspond to the non-relativistic vertex oc [cr x fc]. 
Coupling constants can be, e.g., taken from the analysis 
of the Jiilich model of the hyperon-nucleon interaction 
via the meson exchange [451: Ck*na = , Ck*ny. = 

, and Ck^ns* = — ■ These values account for the 
form-factors used in [45]. In particular, the form-factor 
related to a very soft energy range is responsible for a 
strong suppression of the Ck*ny, vertex. 

Thus, the Landau-Migdal correlation parameters (72) 
can be cast as 


Cof'A 

Co5s 

Co5e* 


'-'KNA 
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if ATS 


KNA CkNS 


3 Catate* 


■ 4 ^’“’''^^( 0 ) + i?AA4P’( 0 ) 

4’™’''=)(0) + i?EE4.’“’'‘’^(0) 

4’’™’''^(0) + i?E*E* 42’“’''^ 


where the first term was introduced in (68), 4’ (0) ~ 

0.24, the second one is equal to = 

|m^/(m^ -b m^.) ~ 0.28, cf. [57], Rhh' = 

Ck*nh Ck*nh'/{Cknh Cknh') with i?AA — 3.7, 

i?EE — 0.04, i?AS — 0.39, and i?E*E* — 0.69. Additional 
factor 2 in 4* compared to 4’ j originates from 
the reduction [<72 x fc] [cri x fc] ^ S 12 . Finally, we esti¬ 
mate the following values of the correlation parameters 
in (72) as 

/;~0.9, 5E-0.1, /;s--0.1, 5 e.- l-2.(77) 

Compared to [15] we obtained a smaller value of the cor¬ 
responding parameter Cq/a — 0.6/m^, since we included 
here the form-factors mentioned above [77]. 


D. Energy of the Lowest Branch of the Dispersion 
Relation at fe =0. 


Fig. 9 illustrates how much baryon-baryon correla¬ 
tions affect the terms of the bare polarization operator. 
For case I we show the energy of the lowest branch of 
the dispersion equation at fc = 0 calculated with the 
polarization operator Ilg^^ constructed from ac¬ 

cording to (64) and the polarization operator ns(w) = 



FIG. 9: Solid lines present the energy of the lowest branch 
of the dispersion equation at fc = 0 calculated with the po- 
larization operators Ilg ' and IIs including effects of baryon- 
baryon correlations. Dash-dotted continuations of solid 
curves demonstrate effect of the filling of the hyperon Fermi 
seas on the s-wave polarization operator. Correlation parame¬ 
ters are taken according to Fig. 8 (right panel), (66) and (77). 
Dashed lines show solutions and Ilg^^ without correla¬ 

tion effects, as in Fig 6, plane (b). The dotted line depicts 
the electron chemical potential. 


4^^(w) -bnP°i®(w, 0), where nP°i®(w, 0) follows from (73) 
with parameters (77). 

At p > pc,H we have to include correlations in the term 
dn^yp in (61). The pole term 0) is included 

in (73) with the help of the replacement (57). The other 
terms <5n4yp’°^ corrected in the 

same manner as the 11® term. However these terms are 
rather small as it is demonstrated by Fig. 7. Therefore 
we omit correlations in them in further. 

From Fig. 9 we see that baryon-baryon correlations 
suppress the s-wave part of the polarization operator, in 
agreement with the statements [53, 54]. It results in an 
increase of the critical density of the s-wave condensa¬ 
tion, from 2.7pq to 4.3po for case I chosen by us for an 
illustration. 


VIII. K- CONDENSATION IN NEUTRON 
STARS 

In sections IV-VII we have constructed the K~ po¬ 
larization operator. Now we use it to study a possible 
instability of the system with respect to a phase transi¬ 
tion into a state with K~ condensate. 

For this aim we first investigate solutions of the K~ 
dispersion relation 

uj'^ — — m\ — Ren‘°*(a;, fc) = 0, 


(78) 
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where the complete polarization operator is given by 

= ns(w)+n;?®(a;,fe)+nP°‘®(a;,fe) 

+ K™H+^Cyp°M- (79) 

It contains the s-wave part, regular p-wave and pole parts 
of the polarization operator given by (64,65) and (73), 
respectively, and the terms determined by the hyperon 
populations (59) and (60). The correlation parameters 
are taken according to (66,67) and (77). 

There are two different possibilities. The K~ con¬ 
densation may occur in the neutron star matter via a 
second-order phase transition or a first-order phase tran¬ 
sition. The dynamics of both phase transitions are quite 
different. Therefore, both possibilities might be realized 
at different physical conditions related to the different 
stages of neutron star evolution. 

In the case of a second-order phase transition, at the 
moment, when the density in the neutron star center 
achieves the critical density the reactions (I) come 
into the game. At this second order phase transition 
the isospin composition and the density may change 
only soothly. For typical times r oc Treact \/^/ \/p — P^c 
the system creates an energetically favorable condensate 
state, Treact is the typical time of weak processes (I). The 
condensate appears within the region where p > p]}. If it 
happens during the supernova explosion, the typical size 
of the condensate region might become of the order of 
the neutron star radius. Due to the energy conservation 
the positive energy is released in such a transition. When 
the condensate region is heated up to the temperatures 
T > Topac ~ (1 2) MeV neutrinos are trapped. At this 

stage the cooling time is determined by the neutrino heat 
transport from the condensate interior of the neutron star 
to its exterior [24]. When the star cools down to smaller 
temperatures, T < Topac, the neutron star begins to be 
transparent for neutrinos. They can be directly radiated 
away. A part of the energy is radiated by photons from 
the star surface. In binary long lived systems the critical 
density in the neutron star center can be achieved by the 
accretion process. Then the transition is characterized 
by the typical large time of the accretion. 

In the case of a first-order phase transition the final 
state might significantly differ from the initial one by its 
isospin composition and the density. Thus this new state 
can’t be prepared in microscopic processes. Too small 
droplets of the new phase are not energetically favorable 
due to a positive surface energy contribution. When the 
density in the star center begins to exceed the value p\ < 
p]} the system arrives at a metastable state. When in 
a fluctuation a droplet of the new phase (with density 
Pc” > Pc) of rather large (overcritical) size is prepared it 
starts to grow. At zero temperature, the probability of 
the creation of such a droplet via quantum fluctuations 
is very small but increases greatly with the temperature 
[24]. Thus, the first-order phase transition occurs most 
likely at an initial stage of the neutron star formation 
or cooling when the temperature is rather large. If the 


density in the center of a star exceeds the value a 
second-order phase transition may also occur. In binary 
stellar systems, where the neutron star slowly accretes 
the mass from the star-companion and the temperature 
is small, the second-order phase transition might be a 
more probable one (depending on the accretion rate). 


A. II Order Phase Transition to the s-wave 
Condensate State 

Let us first analyze possibility of the s-wave K~ con¬ 
densation. 

In Fig. 10, summarizing the results of sections V C, 
VIC, VIID, we show the energy of the lowest K~ branch 
of the kaon dispersion relation (78) for fc = 0 together 
with the electron chemical potential. For the given pa¬ 
rameter choice (see different panels), the crossing points 
of the lines indicate the critical density of the s-wave con¬ 
densation via the second order phase transition. 

We see that for all the models the neutron star mat¬ 
ter is unstable with respect to reactions (1) in the density 
interval (3-^5.2)po in dependence on the choice of the cor¬ 
relation parameters and the parameters of the hyperon- 
nucleon interactions. Recall, that in the framework of our 
model the attraction due to the s-wave kaon-nucleon in¬ 
teraction corresponds to the effective kaon-nucleon sigma 
term 'Skn equal to 150 MeV, which is much smaller than 
the value T,kn — (300 400) MeV used in the previous 

works studying the s-wave K~ condensation. Additional 
suppression comes from the nucleon-nucleon correlations. 
Nevertheless, we found also an extra attraction associ¬ 
ated with (Jn^s, nP°*® and (511;, terms. Thus, until corre¬ 
lations are not included, we support a conjecture of previ¬ 
ous works on the possibility of the s-wave condensation at 
T'c^S ~ ^Po- However, the baryon-baryon correlations ad¬ 
ditionally shift the condensation critical density to larger 
values than those ones discussed in Refs. [5, 7-12]. 

Note that here we checked the necessary condition of 
the s-wave condensation but we did not yet minimize the 
energy of the lowest branch of the dispersion equation, 
u;™“(fc), over |fe|. Thus we still can’t conclude whether 
we deal with the s-wave or the p-wave condensation. 


B. II Order Phase Transition to the p-wave 
Condensate State 

In this section we are going to study the principal pos¬ 
sibility of the p-wave K~ condensation in neutron star 
matter at the assumption of the second order phase tran¬ 
sition. We would like to investigate whether the p-wave 
condensation second order phase transition can occur at 
densities smaller than that for the s-wave condensation, 
and how much such a situation is sensitive to the param¬ 
eter choices for hyperon-nucleon interactions and corre¬ 
lations. 
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FIG. 10: The energy of the lowest K~ branch of the dispersion relation (78) for fe = 0 is depicted by solid lines. The dotted 
lines show the electron chemical potentials. The dashed lines are calculated without baryon-baryon correlations. Different 
panels correspond to the interaction given by the cases I-IV. 


Let us first assume that the s-wave K~ condensation is 
indeed possible at some critical density Pc,s- The lowest 
energy branch of the K~ spectrum at small momenta is 
given by 

w WS + q;(ws) ^s(ws) , (80) 


where ws is, as before, w(fe = 0) for the lowest energy 
branch of the dispersion law given by the solution of the 
equation Wg = + Rens(ws) ^ = 0), and 


~ 1 2 It’s — 

a{uj) 


i9Rens(w, fc) 


= 1 + 


doj 

9Renp(a;, k) 


fe=0 


> 0 , 


dk‘ 


k=0 


The regular part follows from (41), Q!reg(w) = 
bp(oj) Pp/po + bn{uj) pn/po, with the coefficients defined 
in (42,43). Suppression of the regular part Oreg due to 
the baryon-baryon correlations can be taken into account 
according to (65) 


/ \ bp(uj^ bni^UJ^ ‘Zbpi^Uj') brfi{uj')^np (^) 

reg j — 7 ^\ 5 

1 - hp(ljj) bniu!){^np)^iuj) 

y ,,,, _ bp{uj)pp/po 

— /p'. , 

1 - bp{ui) pp Qp'(uj)/po 

Kiu;) = - br.i^)pJPo -, (83) 

1 bn{uj) Pn ^nn (^)//^0 


If a(a;s) < 0 at pc,S; then instead of the s-wave conden¬ 
sation we, actually, have the p-wave condensation at a 
somewhat smaller density. The aim of this section is to 
find the value Q;(a;s) in (80) at the critical point of the 
s-wave condensation i.e., when us = Pe- 

Taking into account the p-wave kaon-baryon interac¬ 
tion which we have determined in Sect. V, we find 


a{u) = 1-1- apole(w) -I- Q:reg(w) . (81) 


Without baryon-baryon correlations, the contribution of 
the pole part is 


( 0 ) 


pole 


dk- 
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From (39,40) we have 


— ^KNH (vnH + VhN ^ 

= C'lrATE* ^^E*‘J’iE* (w, 0), 7L = {A, S}, . 

The term (j)^j^{—u) accounts for the contribution of 
the hyperon Fermi sea. Once the baryon-baryon corre¬ 
lations are included in the pole part of the polarization 
operator according to (73), ctpo^e is to be replaced by 

apoie(a;) = T|2l^enP°“^(u;,fe) . 

ak fc=o 


/p'\ 

Although the ^7, (w)-functions can be evaluated as in 
(67), we will treat them here as free energy-independent 
parameters ^u ', in order to investigate the sensitivity of 
the results to them. 

In Figs. 11-12 we show values of Q!poie(lie) (solid lines) 
and —1 — Oreg (dash lines), calculated as a function of 
Pc,s for various baryon-baryon correlation parameters 
and for cases I-IV determining the hyperon couplings 
in (2). The S hyperon contribution is proved to be very 
small and most part of the strength is due to the Kp~^ 
and E*n“^ contributions. Filling of the hyperon Fermi 
seas is not incorporated (we artificially suppress terms 
oc in nP°^®). As we mentioned, embedding S* into 
the mean-field model (2) is quite uncertain due to the 
absence of any empirical constrain on the coupling con¬ 
stants. Therefore we consider two different cases. In 
the upper plot in Fig. 11 we assume that E* couples to 
the mean-field with the same strength as the E-hyperon 
(Vs* = Ve, ^ctE* = Scte), whereas in the lower plot 
in Fig. 11 we detach E* from the mean field potentials 
(Ve* = 0, 5crE* = 0). 

At the crossing point of the solid line with the corre¬ 
sponding dash line we have Opoie = — 1 — Oreg and there¬ 
fore a = 0, that means that given density is the criti¬ 
cal density for both the s- and the p-wave condensations. 
When at the given value of pc,s the solid line is below the 
corresponding dash line we have Opoie < — 1 — Oreg. This 
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FIG. 11: The quantities apoie(lie), tus = /Xe, (solid lines) and —1 — Oreg (dash lines) are shown vs. critical density of a s-wave 
kaon condensation for various correlation parameter choices. Labels show values of /)( = Solid lines are drown for /)(j; = 0, 
and for /a = 5s: = 0,0.5, 1. Dashed lines are drown for three choices of the correlation parameters used in calculation of areg, 
Cnn = Cpp = ^np = 0,0.5,1 (from the upper line to the lower one). Cases (I-IV) correspond to those in Fig. 1. In the upper 
plot the E* contributions are taken with = gs, Ve* = Vs and 5 o-e* = So-e. In the lower plot the E* is detached from 
the mean-field potentials Ve* = 0 and 5ctE* = 0. Hyperon populations are included in mean fields but not in the polarization 
operator. 


means that a < 0 and for such a parameter set K~ con¬ 
densation occurs in the p-wave state at somewhat smaller 
density than the assumed s-wave condensation. 

The contributions from the hyperon Fermi seas are in¬ 
cluded in Figs. 12. From Figs. 11 and 12 we observe 
that S* contributes essentially to the K~ polarization 
operator increasing attraction in the p-wave. The most 
favorable case for the p-wave condensation is realized if 
S* is detached from the mean field potentials. Compar¬ 
ing Figs. 11 and 12 we also see that the effects from the 
hyperon population on the p-wave terms are large. 

In Fig. 13 we show the results of our calculations for 
the values of the correlation parameters, which we have 
evaluated in Sect. VII. For the E* detached from the 
mean-field potentials our model predicts a preference of 
the p-wave condensation. The p-wave condensation may 
occur already at p < 3po for the cases II, IV, and for 
p ~ 4po and p ~ 4.5po for cases III and I provided the 
s-wave softening of the spectrum is also rather high. For 
S* coupled with the same strength as E the s-wave con¬ 
densation might be preferable compared to the p-wave 
condensation. Since a variation of other parameters, be¬ 
sides E* strength, is certainly also allowed the results are 


more uncertain in the latter case, cf. Figs 11 and 12. 


C. I Order Phase Transition to the K Condensate 

State 

In the previous two sub-sections we have studied the 
possibility of the K~ condensation assuming that it oc¬ 
curs by a second-order phase transition. Now we will 
investigate properties of K~ excitations in the baryonic 
matter of different particle compositions, to hgure out 
whether an abrupt ( of the first order) phase transition 
into a state with a new particle composition and another 
baryon density is energetically favorable. 

We shall consider: (i) nucleon-hyperon matter (NHM) 
composition which we have discussed above, see Fig. 1; 
(ii) neutron-proton matter (NPM) composition, the mat¬ 
ter consisting only of protons and neutrons assumed to be 
in /3-equilibrium; (iii) isospin-symmetrical nuclear matter 
(ISM) composition consisting of protons and neutrons 
with pp = pn and leptons, compensating the electric 
charge of the protons; and (iv) the proton-enriched mat¬ 
ter (PEM) consisting of protons and neutrons with the 
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FIG. 13: Values Opoie (solid lines) and —1 — areg (dash lines), cjs = Ate, calculated for the correlation parameters (77). Dotted 
lines present calculations with E* detached from mean held potentials. Planes correspond to the cases I-IV of hyperon couplings. 
Possibility of the hyperon population is also incorporated in the polarization operator. 


concentration Yp = Pp!p = 0.7 and the charge compen¬ 
sated by the leptons. Thus in cases (ii)-(iv) we switch 
off the hyperons in our mean-field model (2) and in cases 
(iii), (iv) we also freeze the value of the proton concen¬ 
tration. For each case we shall calculate the total energy 
of the system with and without K~ condensation. The 
idea behind this procedure is as follows. We consider 
the given configuration with the corresponding proton 
concentration and the condensate field as the most ener¬ 
getically favorable configuration until we did not prove 
the opposite. The full reaction balance is assumed to be 


recovered, thereby. Below we show that the first order 
phase transition occurs from NHM starting at the density 
p\ < P^c to approximately ISM at p\^ > p\. 

In Fig. 14 we show the energy of the lowest K~ branch 
of the dispersion equation (78) Yp, p) minimized 

with respect to the momentum, as a function of the 
baryon density for different proton concentrations. We 
see that the more protons exist in the matter the smaller 
is the value of the density at which the minimal kaon 
energy w™'"(fem, 1^, p) meets the electron chemical po¬ 
tential. Therefore if the energy gain due to the conden- 
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FIG. 14: The energy of the lowest K~ branch Y^, p) minimized in |fel, and the electron chemical potential, calculated 

for the NHM, the NPM, ISM and PEM. In the right and left panels the results are presented with and without baryon-baryon 
correlation effects, respectively. In the NHM case interaction forces of case I are used and the E* contributions are taken with 
Se* = Vs* = He and = 5<tE. 


sation were large enough to compensate an energy loss 
in the fermion kinetic energies the system would undergo 
a first-order phase transition from the NHM state to the 
state with a proton enriched composition and a different 
density. 

Let us first compare the energies of NHM, NPM, 
ISM and PEM taking into account a possibility of the 
K~ condensation in each case. For densities, when 
p) > pe, there is no condensation and the 
energy density of the system is given by (4). When, at 
given Yp and p, we have uj’^'’^{kjn,Yp, p) < pe, the K~ 
excitations appear, partially replacing the leptons. These 
excitations occupy single state and form a K~ conden¬ 
sate. The electron chemical potential is fixed now as 
Pe = Yp,p). The energy density of the system 

with the K~ condensate reads: 

B 

+ EiiL0’^^’^{krn,Yp,p))+El^J^, (84) 

l—fi e 

where E^^^^ is the energy density of the condensate field 
related to the mean field Lagrangian 

Lk = {io^-IV{ijj,k))\(j)k\’^. (85) 

The (j)k is the K~ condensate mean field with the wave 
vector k. This field component should be found from 
the minimization of the appropriate thermodynamical 
potential. For simplicity we use here the linear model 
neglecting the higher order terms, as oc \(j)k 1^, which 
incorporate an effective kaon-kaon interaction in dense 
baryonic matter. The effective kaon-kaon interaction de¬ 
pends on the structure of the mean field. In absence of 
the non-linear effective kaon-kaon interaction, within the 
variational procedure we obtain the field of the running 


plane wave type, ipk = exp(—p)i + ifemi'). 
Using that the dispersion relation (78) is fulfilled for 
k = km and ui = u!™^’^{km,Yp, p), and that the den¬ 
sity of the charged kaon condensate px is fixed by the 
electro-neutrality condition 

PK = QB PB - Pe - P^- , 

B 

we find the energy density of the kaon condensate equal 
to 

Yp, p){^^qBPB- Pe- P;.-) ■ 

B 

In Fig. 15 we show the energy per baryon in the vari¬ 
ous baryonic systems, NHM, NPM, ISM, and PEM {Yp ~ 
0.7), with and without the kaon condensate (dashed and 
solid lines, respectively). We see that for the density 
Pc — 3 po without baryon-baryon correlations and ~ 4 po 
including correlations, the condensate state of ISM be¬ 
comes to be energetically favorable compared to NHM. 
The transition to the new more symmetric isospin con¬ 
figuration would increase the Fermi energies of the lep¬ 
tons (the latter ones are needed to compensate a larger 
charge of protons), but reduce the symmetry energy of 
the nuclear matter and the total Fermi energy of nu¬ 
cleons. Without the K~ condensate, the energy loss is 
larger than the gain and the system chooses the composi¬ 
tion shown in Fig. 1 for the given EoS, cf. that in Fig. 14 
the solid lines labeled ISM are lying above the NPM lines. 
With account for a possibility of the K~ condensation 
one additionally gains on the kinetic energy of the lep¬ 
tons, which are replaced by kaons. The latter energy gain 
is large enough, leading to the preference of the ISM. The 
hyperon Fermi seas are not filled. It is more energetically 
preferable to create extra K~ condensate mesons instead 
of a filling of a Fermi sea of the hyperons. From Fig. 15 
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FIG. 15: Total energy per baryon of the matter with different particle compositions: without a K~ condensate (solid lines) and 
with the condensate (dashed lines). Abbreviations are the same, as in Fig. 14. The curves in the right and the left panels are 
drawn with and without inclusion of correlations. In the NHM case interaction forces of case I are used and the E* contributions 
are taken with — gs, Vs* = 14; and g^s* = ga-s- 


we see that PEM {Yp ~ 0.7) has a larger energy than 
ISM. As it is also seen from Fig. 15 the resulting isospin 
composition can be only slightly above Yp = 1/2. Fur¬ 
ther on we neglect this difference considering ISM as the 
final configuration. 

In Fig. 16 we plot the lowest branch of the K~ excita¬ 
tion spectrum Wmin(|fc|) in the ISM at various densities. 
In the left panel baryon-baryon correlations are switched 
off and in right panel, switched on. We see that for rather 
large densities (> (2.5-h 3.5)po) the spectra have minima 
(dots in Fig. 16) at finite values of the kaon momentum 
km 7 ^ 0. It signals that the transition from NHM to 
a dense ISM would occur as the first-order phase transi¬ 
tion into the state with the p-wave kaon condensate. The 
calculations shown in Figs. 14-16 are done for = 9's, 
Vs* = 14 and ga-s* = ffo-s- The results obtained with the 
E* detached from the mean field potentials are checked 
to be of minimal difference. 

In the assumption that the surface tension is larger 
than a critical value the initial and final state densities 
are determined by the double-tangent Maxwell construc¬ 
tion [58, 59]. In Fig. 17 we show the double-tangent con¬ 
struction between NHM and ISM states. We see that the 
critical density for the beginning of the first-order phase 
transition is equal to p], ~ l.dpo without correlations and 
~ 2.1po with correlations. The critical densities of the fi¬ 
nal state are p\^ ~ bpo and ~ bpo, accordingly [78]. 

If surface tension is smaller than the critical value, the 
phase transition results in the mixed phase [58-61]. Then 
the local charge-neutrality condition is relaxed being re¬ 
placed to the global charge-neutrality condition. In this 
case the critical density for the appearance of the kaon 
condensate droplets within mixed phase is still smaller 
than that value given by the Maxwell construction. The 
presence of the mixed phase may have interesting observ¬ 
able consequences, see [62] and references therein. 


Thus, relying on the analysis above, we argue that the 
actual critical density of the first-order phase transition 
can be even smaller than 2pQ and that this transition 
occurs into the p-wave condensate state. 


IX. CONCLUSION 

In this work we constructed the K~ polarization oper¬ 
ator in dense baryonic matter of arbitrary isotopic com¬ 
position, including both the s- and p-wave -baryon 
interactions, and, using the relativistic mean field model 
to describe the baryon properties. We applied the de¬ 
rived polarization operator to the issue of the s- and p- 
wave K~ condensations in neutron star interiors. We 
considered two different models of the equation of state, 
cf. Sect. H and Appendix A below. Finite temperature 
effects can be easily incorporated in our general scheme. 

To describe the kaon - nucleon interaction we used the 
kaon-nucleon scattering amplitude obtained as a solu¬ 
tion of the coupled-channel Bethe-Salpeter equation with 
the interaction kernel derived from the relativistic chiral 
SU(3) Lagrangian with the large N,. constrains of QCD. 
We calculated explicitly the pole terms of the K~ po¬ 
larization operator related to exci¬ 

tations with K~ quantum numbers and analyzed effects 
of the filling of the hyperon H — (A, E, 5) Fermi seas 
at densities above the hyperonization point p > pc,H — 
(2.5t3)po- 

In Fig. 6 we compared the s-wave regular part of our 
polarization operator with the simplified form, which is 
widely used in the literature. The E-term extracted from 
this comparison (E ~ 150 MeV) is found to be two-three 
times smaller compared to that allows for the s-wave K~ 
condensation in ordinary neutron star matter composed 
mostly of neutrons. However, we found the essential at- 
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FIG. 16: The energy of the lowest branch of the dispersion equation (78) as function of the momentum with the polarization 
operator (79) for the ISM at various densities (labels on the curves show the densities in po). Calculations, done without and 
with baryon-baryon correlations, are shown on left and right panels, respectively. Full circles mark the position of the minimum 
on the spectrum. 




FIG. 17: Energy densities of NHM and ISM with and without the K~ condensate shown by dashed and solid lines, respectively. 
For NHM the forces of case I are used. Calculations with and without baryon-baryon correlations are shown in left and right 
panel respectively. The dash-dotted lines represent the double-tangent Maxwell constructions between two phases. 


tractive support from the hyperon exchange terms of the 
p-wave scattering amplitude contributing to the s-wave 
part of the polarization operator, see (48). Inclusion 
of these terms, which were omitted in previous works, 
makes a second-order phase transition to the s-wave K~ 
condensate state possible (at densities ~ Spo in given 
model, when the correlation effects are not included yet). 

We evaluated baryon-baryon correlation parameters 
and corrected all the s- and p-wave terms of the polariza¬ 
tion operator, accordingly. Inclusion of the correlations 
pushes the critical point of the second order phase tran¬ 
sition to the s-wave K~ condensate state up to larger 
densities (4 d- 5)poj cf. Figs. 10 and 20). We studied 
how much the results are sensitive to a variation of the 
correlation parameters. We estimated (see Appendix D) 
feed-back quantum fluctuation effects arguing that their 
contribution is not too large at the low kaon energies un¬ 
der consideration, and at zero temperature to the first 
approximation they can be neglected. 


Our next observation (see Appendix D) is that at fe = 0 
the imaginary part of the pole term of the polarization 
operator is finite only in a rather narrow interval of kaon 
energies. Would the electron chemical potential cross 
the K~ branch within this energy interval, the s-wave 
condensation would not occur. However this possibility 
is not realized for the parameter choice used in our model. 

We checked the possibility of a second-order phase 
transition to the p-wave K~ condensate state. We 
showed that, in the vicinity of the critical point of the 
s-wave K~ condensation, the p-wave part of the polar¬ 
ization operator, induced mainly by A-proton holes and 
S*-nucleon holes and some regular terms, is rather large 
and attractive. It may change the sign of the momen¬ 
tum derivative of the energy at the lowest K~ spectrum 
branch at origin. If occurred, it would mean that we have 
p-wave condensate instead of the s-wave one at somewhat 
smaller density. We demonstrated that this statement, 
although being rather model dependent, holds for a wide 
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range of varying parameters. The results essentially de¬ 
pend on the unknown value of the strength of the S* 
mean field potential. In most favorable case, when E* is 
detached from the mean field potentials, the second or¬ 
der phase transition to the p-wave condensate state may 
occur even at p ~ 3po (with correlations included), cf. 
Fig. 13 and 21. This result is also sensitive to the details 
of the equation of state and the parameterization of the 
hyperon-nucleon interaction. For the equation of state 
with parameters (Al) the critical density is increased 
compared to that one calculated with parameters (5). 

We discussed the possibility of a first-order phase tran¬ 
sition to a K~ condensate state. We calculated the en¬ 
ergies of the baryonic matter with different compositions 
with and without the inclusion of the K~ condensation. 
We found that with account for a K~ condensation, the 
isospin-symmetrical neutron-proton matter is more en¬ 
ergetically favorable than the standard nucleon-hyperon- 
lepton matter at densities 3 4po (depending on the 
values of parameters of baryon-baryon correlations). The 
hyperon Fermi seas are melting at this phase transition. 
Hyperons are replaced by nucleons and electrons are re¬ 
placed by the condensate K~ mesons. We demonstrated 
that in dense isospin symmetrical np-matter the K~ ex¬ 
citations are condensed in the p-wave state. With the 
help of the Maxwell construction we found that the crit¬ 
ical density of the beginning of the phase transition is 
about 2po with the baryon-baryon correlations included, 
cf. Figs. 17 and 24. The final state density is about 
(5 6)po- Appearance of such a strong first order phase 

transition may have interesting observable consequences 
as blowing off a part of the exterior of the neutron star, a 
strong neutrino pulls, a gravitational wave, a strong pul¬ 
sar glitch, etc. These consequences have been previously 
discussed in relation to the first order phase transition 
to the pion condensate state [24]. Here we may expect a 
stronger energy release compared to the pion condensate 
phase transition since typical energy scale rriK ^ TOtt- 

Our derivations can be helpful not only for description 
of neutron star interiors, but also for discussion of kaonic 
effects in other nuclear systems, as atomic nuclei and 
heavy ion collisions. Therefore, of particular interest is 
the further more detailed analysis of the p-wave effects 
on K~ spectra in nucleus-nucleus collisions motivated by 
present SIS and SPS experiments and the future SIS200 
program. 

In spite of a number of new effects was incorporated 
in our scheme, some other effects might be also impor¬ 
tant. Present calculations still suffer of many uncertain¬ 
ties, most of which are due to the lack of experimen¬ 
tal information, e.g. on the coupling constants, the ab¬ 
sence of unambiguous way for going off-mass shell and 
the lack of study of more complicated in-medium fluc¬ 
tuation effects, which we just roughly estimated in the 
present work. Among them, there are the pion softening 
effects [24], which can essentially affect the results at fi¬ 
nite temperatures, and the contribution of the non-linear 
meson-meson interaction. The latter may partially sup¬ 


press the condensate contribution to the energy beyond 
the critical point. 

Diminishing of the uncertainties needs further theoret¬ 
ical and experimental work. 
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APPENDIX A: VARIATION OF THE 
PARAMETERS OF THE BARYON 
INTERACTION 


In this section we investigate how much the results of 
Sec. VHI are sensitive to the particular choice of the EoS. 
For comparison with the results obtained within the rela¬ 
tivistic mean-field model with the parameters (5) we use 
here the EoS in parameterization [63], which is a good 
fit to the optimal EoS of the Urbana-Argonne group [26] 
up to a 4 times nuclear saturation density, smoothly in¬ 
corporating the causality limit at higher densities. 

The parameters of the mean-field model are adjusted 
to the following bulk parameters of the nuclear matter 
at saturation: po = 0.16 fm“^, binding energy Abind = 
— 15.8 MeV, compression modulus K =250 MeV, sym¬ 
metry energy Ogym = 28 MeV, and the effective nucleon 
mass m*j^{po) = 0.8 niN- The corresponding coupling 
constants of the Lagrangian (2) are then as follows: 


9wN 


= 91.2506, 


2 2 
9^n rriM 


= 195.597, 


2 2 
9pN rriN 


= 77.4993, 6 = 0.0867497, 


c= 0.0805981. 


(Al) 


In Fig. 18 we show energies for the nucleon isospin- 
symmetrical matter (ISM) and the pure neutron matter 
for two choices of the mean-field EoS, the model (5), sim¬ 
ulating a softer EoS and (Al) simulating a stiffer Urbana- 
Argone EoS. In spite of the parameters (5) and (Al) es¬ 
sentially deviate from each other, the energies and other 
thermodynamic characteristics of the neutron star mat¬ 
ter are rather closed to each other in both the parameter 
choices in the absence of a K~ condensate. For the ISM 
case the EoS with parameters (Al) is significantly stiffer 
than the one calculated with parameters (5) at p U 3po. 

In Eig. 19 we show concentrations of the baryon species 
in neutron star matter corresponding to EoS given by the 
choice (Al) for four choices of the hyperon-nucleon inter¬ 
action (cases I - IV) which we are using in this paper. 
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FIG. 18: EoS for ISM {N = Z) and neutron {Z = 0) matter 
calculated with the parameter set (5) (dotted lines) and with 
the parameter set (Al) (dashed lines). The solid lines show 
the EoS from Ref. [63]. 


We see that the critical density of the hyperonization is 
approximately 3po for all the choices. General trends are 
the same as the ones shown by the corresponding curves 
discussed in the main text, see Fig. 1. The most essen¬ 
tial difference is that the proton concentrations given by 
Fig. 19 are smaller than those in Fig. 1. This should have 
consequences for the neutrino cooling of the neutron star. 
Indeed, when the proton concentration enlarges ll-hl4% 
the efficient cooling mechanism is switched on due to al¬ 
lowance of the direct Urea processes p —> n + e + v. In 
principle, this difference may allow to select a more ap¬ 
propriate EoS in the future. 

Fig. 20 demonstrates the energy of the lowest K~ 
branch of the dispersion equation (78) for k = 0 calcu¬ 
lated with the mean field model with parameters (Al). 
Dashed lines are computed without inclusion of correla¬ 
tions and solid lines, with inclusion of correlations. Com¬ 
paring this figure with Fig. 10 presenting the same cal¬ 
culation but with the parameter choice (5) we see that 
the critical densities are increased for all the cases I-IV 
for the model (Al). The second order phase transition to 
the s-wave K~ condensate state occurs within the den¬ 
sity interval (3.5 6)po depending on the choice of the 
parameters. For the model (5) the corresponding interval 
of critical densities was (3 -h 5.5)po- 

In Fig. 21 we demonstrate the possibility of the sec¬ 
ond order phase transition to the p-wave state within 
the model with parameters (Al) at the assumption that 
the system is in the vicinity of the critical point of the 
s-wave condensation. Solid curves relate to the model 
when E* is attached with the mean field potentials with 
the same strength as E and dashed curves, when E* is 
detached from the mean field potentials. Comparing the 
curves with those curves presented in Fig. 13 for the pa¬ 


rameter choice (5) we see that the general trends are the 
same in both cases. The second order phase transition 
to the p-wave state may occur in the model (Al) for 
P > (4.5 -h 5.5)/9o in cases II - IV and it does not occur 
for case I up to 6po. 

Figs. 22, 23, 24 study the possibility of the first- 
order phase transition for the EoS with parameters (Al). 
Again we see the same trends as given by corresponding 
Figs. 14, 15, 17. From Fig. 22 we see that the energy 
of the PEM crosses the electron chemical potential at a 
smaller density compared to ISM and NHM cases. Fig. 23 
demonstrates that the energy of the ISM with the con¬ 
densate (with correlations included) becomes to be lower 
than the energy of the NHM for p > ApQ. This value (dpo) 
is about the same as that for the EoS with parameters 
(5). 

In Fig. 24 we show the double-tangent construction 
for the EoS with parameters (Al). The first order phase 
transition starts at the density p], ~ 2.5po (with corre¬ 
lations included). This value is only slightly larger than 
that (2.1po) given by the EoS with parameters (5). The 
final state density is p\^ ~ 5.5/9o, he. slightly less than 
the value (5.9/9o) given by the model (5). 

Thus the main conclusion, what we can do from the 
above comparison of the results obtained for two EoS 
with parameters (5) and (Al), is that all the general 
trends in the behavior of the kaon condensation are the 
same. The critical densities of the s- and p-wave conden¬ 
sations are only slightly higher in the framework of the 
model with parameters (Al). 


APPENDIX B: NON-EQUILIBRIUM GREEN’S 
FUNCTIONS 


Two-point (2P) functions (a Green’s function or a 
self-energy) are introduced within the Schwinger-Baym- 
Kadanoff-Keldysh approach as, cf. [64], 


iF{x,y) 


iF {x,y) iF +{x,y)\ 
iF+-{x,y) iF++{x,y) J 


' lrA[x)B{y)) 

^ ll{x)B{y'^ 


^lB{y)A{x^ ^ 
(T-^A{x)B{y)) ^ 


(Bl) 


where T and T~^ are the usual time and anti-time order¬ 
ing operators. Note that, compared to Green’s functions, 
“± =f” self-energies would have extra sign — since they 
include vertices V~ = —iVo and U+ = +iVo. 

Not all four components of F are independent. The 
following relations between non-equilibrium and usual re- 
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FIG. 19: Concentration of baryon species in neutron star matter for EoS (Al). Panels (I-IV) correspond to four choices of the 
hyperon-nucleon interaction parameters specified in text. 



FIG. 20: The same as in Fig. 10 but for the parameters (Al). 


tarded and advanced functions are fulfilled 


functions in the Wigner representation, 


F^{x,y) = F—{x,y) - F-+{x,y) 

= F+~{x,y) - F''~+{x,y) 

:= 9{xo - yo) {F+-{x,y) - F-+{x,y)) , 
F^{x,y) = F—{x,y) - F+-{x,y) 

= F~+{x,y) - F''~+{x,y) 

:= 9{yo - xo) (F"+(x,j/) - (x, j/))(p2) 


= J ^e-FP<)FiX + ^,X-0 


(B4) 


where 9{xo—yo) is the step function of the time difference. 
For such 2P-functions complex conjugation implies 


[iF '^{x,y))‘" = iF ~'~{y,x) => iF {X,p) = real, 
{iF^ {x,y))‘" = iF^~{y,x) iF^ (X,p) = real, 

{iF {x,y))‘^ = iF++{y,x) 

^ {iF—{X,p))*=iF++{X,p), 

{F^{x,y)f = F^{y,x) 

=> {F^{X,p))* = F^{X,p), (B3) 


where the right parts specify the properties of the 2P- 


with X = i(x + y) and p = (e,p). The conjugation 
operation is defined as a (...)*^ = (...)* for bosons and 
(...)^ = 7o(...)^7o for fermions. 


We denote the fermionic Green’s function and self¬ 
energy as G“’^ and respectively, and the bosonic 

ones as £)“’*' and 11“’*'. The hats on the fermionic 2P- 
functions point on their spin structure, e.g., — 

[/i — TOf — 


For systems in equilibrium, the “—h” Green’s func¬ 
tions, the spectral functions and occupation numbers are 











E Jp-m^ [Mev] co'^'"(k^) [MeV] 
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FIG. 21: The same as in Fig. 13 but for the parameters (Al). 




FIG. 22: The same as in Fig. 14 but for the parameter set (Al) 



FIG. 23: The same as in Fig. 15 but for the parameters (Al). 
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FIG. 24: The same as in Fig. 17 but for the parameters (Al). 


related by the Kubo-Schwinger-Martin condition: 

iG~+{p) = -Af{p)n\e ), 

tG+"(p) = Ai{p) [1 - n\e )], 

tD"+(p) = Ah{p)n^{e), 

iD+-{p) = A^a{p)[l + n^{e)\, (B5) 

fE"+(p) = 7 f(p) n‘'(e), 

fE+"(p) = -%{p) [1 - n^(e)], 

m-+{p) = -7b(p) n^{t) , 

tn+"(p) =- 7 b(p) [1+ n’"(e)], (B6) 


functions as Gf{p) = + rm) Gf{p) and Af(jp) = (ji + 

m{)A[{p). In the quasiparticle approximation ( 7 f ^ 0) 
we have 


Ai{p) 


ReGf(e,p) 


2TTS[e'^ - elip) - ReE^(e,p)] 
Zl-Sie-Cj,), 

1 4 

2 6p 6 6p 


l-dReE^/de'^ 



(B9) 


where Ai{p) = — 2ImGf(p), Ah{p) = —21mD^{p) are 
the fermion and the boson spectral functions, %(p) = 
-2ImEf(p), 7 b = -2Imn^ are the corresponding 
widths, and 

= {exp[(e - Pf,h/T] ± 1}“^ (B7) 


are fermion/boson occupation numbers, with /if^b stand¬ 
ing for the fermionic and bosonic chemical potentials. 

In the quasiparticle approximation { 7 b —> 0 in the 
bosonic Green’s functions) we have for bosons 


^b(g) 


ReDl^iq) 


2TTS[qQ -q^ - ml- Ren^(go,g)] 


E 


zY 


90 - 


(B 8 ) 


where Cp obeys the dispersion equation Cp = el{p) + 
ReE^(ep,p) with eo(p) = \/rril -f p^. The self-energy 
E^ includes averaging over the spin structure, E^ = 
i Tr{E« • + mf)} 

The mean-field solutions (3) of the equation of motion 
for a baryon B, following from (2), can be parameterized 
by the self-energy E'^(e, p) = 2 Vb e + — 2mN g<T cr + 

9aB Then we have Cp = Eb{p), with Eb{p) defined 
in (3) and Z^ = [1-Vb/Eb{p)]~^ = l + VB/eB{p) ~ 1. 

The in-medium distributions of particles of given 
species over the momenta are [79] 


d^iVb 

d3Xd3/c/(27r)3 

d^N{ 

d?Xd?k/ (27r)3 


r°° Af 

^2eAb(n^ + 1/2), (BIO) 

I (Bll) 


where Z^ = l/[ 2 qo — (q) are quasipar¬ 

ticle normalization corresponding to the given spectrum 
branch w(,(< 7 ). 

As a step to a non-relativistic limit it is convenient to 
approximate the spin structure of the fermionic Green’s 


The terms t 1/2 are due to quantum fluctuations and 
should be properly renormalized with the corresponding 
subtraction of the vacuum pieces. The contributions cor¬ 
responding to the integration over the negative energies 
are to be included in the distribution of antiparticles. 
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APPENDIX C: IMAGINARY PART OF THE 
LINDHARD’S FUNCTION 


APPENDIX D: EVALUATION OF 
FLUCTUATION TERMS 


The imaginary part of the Lindhard’s function (28) is 
determined by the integral 


i(ppi) 


Im$j_f/(a;,fc) = - 


m 


N 


dE6 




27r|fc| 

- {2ujE + io‘^ 


For brevity of notations we omitted here the vector po¬ 
tentials of baryons which can be recovered at the very end 
by the gauge invariant substitution to ^ to + ViH ■ The 
area where the imaginary part is finite is determined by 
the inequality 

4 {E^ - -{2ujE + uj^ -k"^ > 0, 

which has solutions E > E^fj{uj,k) and E < E~^{uj,k) 
for < k^, and E~^{uj,k) < E < E^^{u),k) for < 
(m^ — and > {m\ + -I- k^ with 


Let us discuss corrections to the kaon polarization op¬ 
erator which we did not take into account yet. We 
mean the fluctuation terms, which describe changes of 
the baryon propagation due to radiation and capture of 
virtual mesons, analogous to the Lamb shift in QED. 
Calculating the corresponding diagrams we shall use the 
Schwinger-Baym-Kadanoff-Keldysh (SBKK) technique, 
albeit, for systems at equilibrium and zero temperature 
the standard Feynman technique is equally appropri¬ 
ate. In the SBKK technique together with usual causal 
Green’s functions G there appear time-disordered ob¬ 
jects G and G^ which have meaning of Wigner’s den¬ 
sities, cf. (Bl). In general non-equilibrium case these 
Wigner’s densities obey the generalized kinetic equa¬ 
tion [64], being therefore very convenient for description 
of the fluctuative effects. Some relations between non¬ 
equilibrium Green’s functions were introduced in Ap¬ 
pendix B. We will estimate fluctuation effects to the 
baryon self-energies and feedback from the kaon modifi¬ 
cation in medium to the polarization operator. 




{-LO E{k‘^) ± jfcj Q{k‘^ 

\/w2 - fe2 


1 . 


Correction of the Baryon Green’s Functions Due 
to Kaon Fluctuations 


One can check that E-^{u),k) < m*j^, if < k^, and 
E~^{uj,k) > if > k^. Therefore, in the case 
a;2 < k^ we obtain 

lm^]^{u,k) = - E+j{uj,k)^ 

X 0(^e,{pF^) - E+^{u;,k)y 

The condition ei{pFi) > E^^{uj,k) imposes the following 
restriction on the frequency: 


eniPFi - |fc|) < w + ei{pFi) < eniPFi + |fc|) • (Cl) 
In the other case -I- k^ we find 

Im$^^(w,A;) = -0(^€^(pF^) - E-fj{uj,ky 

For k^ < ui'^ the condition ei(pFi) < E~^{u),k) puts the 
same restriction on frequency, as (Cl). Within this in¬ 
terval we always have €i{pFi) > E^^{u!, k) and hence 


Im$-2 (w, k) 


X 


^ \/w2 - 

^(e^iPF^) - E-jj{uj,k)y 


In Refs. [19, 20, 39-41, 43] it was argued that the 
fluctuation contributions to the hyperon self-energy are 
essential at low densities and at the kaon energies not far 
away from the mass-shell. This is, particularly, due to the 
presence of the dynamically generated A(1405) resonance 
close the kaon-nucleon threshold. This resonance domi¬ 
nates the s-wave kaon polarization operator at w ~ ttik 
and is very sensitive to the Pauli-blocking effect and to 
the modification of the kaon spectral density. At lower 
kaon frequencies, typical in our case, the influence of the 
A(1405) resonance is small. Thus, there remains to an¬ 
alyze the self-energy contributions of A(1116), 11(1195) 
and S*(1385) resonances. Here we will demonstrate that 
in the low energy region and at the proton-neutron den¬ 
sities of our interest these contributions are rather sup¬ 
pressed. 

Let us first show this on example of the diagram 

K 

This is the self-energy insertion to the full hyperon 
Green’s function due to the KN intermediate states. In 
(Dl) we draw full vertices and the nucleon line represents 
the full Green’s function. To be specihc let us concen¬ 
trate on the kaon-proton self-energy insertion for the A 
hyperon. This corresponds to H = A, K = (accord¬ 
ing to selected arrow direction), N = p in the diagram 
(Dl). 
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Within the SBKK diagram technique, using relations 
(B2) between the retarded and , +” Green’s functions 
and self-energies we obtain, cf. [65], 




V 


G- + (p + q)D^(q) 


+ G^(p + q)n~+(q) - iG- + (p + q)2Imn^(q) 


+ G^(p + q)D^(q) 


V 


d^q 

(27r) 


4 • 


(D2) 


Here —iV = —and —iV^ are the full and bare 
” vertices, and iV = iVq F, iVq are the corresponding 
“-b” vertices, F is a scalar form-factor, which includes 
short-range nucleon-baryon correlations and simulates, 
thereby, the difference between the bare and the full ver¬ 
tices. 

The last term in (D2) vanishes identically, since the 
both retarded Green’s functions have the pole in the same 
complex go-semi-plane. 

To separate the contributions from particles and 
antiparticles K~, we may use the following decomposi¬ 
tions of the retarded Green’s functions and the Wigner’s 
densities 


D^{q) = 9{qo)D^+iq)+ei-qo)Di.i-q), (D3) 
D~+{q) = 0{q^)D-+{q)+0{-qo)D+Z{-q). (D4) 

This allows to reduce integration over go in (D2) to the 
positive values, for A self-energy we have 

E«(p) = iTr{S«(;i + ml)} = / 

X i^iG~~^{p + q)[D^+iq)-2ilmD^+{q)] 

+G^(p + q)^D-tiq) 

+zG-+ ip - q) [D^- {q)-2i Im _ (g)] 
+G’^ip-q)^D-^t{q)], (D5) 


where we used notation 


= V,^V\q), 


Vi = C 


2 

KNA 


X 


1 ( ii> + mX) 

2 \ 2m\ 


^75 


ii>-i\- m%) 


2m 


N 


^75|.(D6) 


we shall see below non-quasiparticle corrections to the 
intermediate fermion Green’s functions, although exist, 
produec small contributions to the self-energies at kaon 
energies of interest. 

We first evaluate ImS^. Using (B5), (B8), (B9) we 
get from (D2) 


ImS^(e,p) f d^g ^.2 2 t .2 

2 ml 4^27 (27r)2 0^ 

’*[l-nP_q](5(A!_(e,p,q)) 
+Z^^’*nP+q(5(A!^(e,p, q))|, 

A±(e,P,9) = e±uj]^±iq) - Epip±q), (D7) 


where i runs over all the K~ and AT'*" branches. Here 
uj^iq) are functions of q^, and = 0(pf,p — \p |) is the 
proton occupation function. Following Appendix B we 
put Z{ = 1. We used also that occupations of real K'^ 
and K~ mesons are absent at T = 0 (we do not consider 
here the processes on the kaon condensate field). 

As it is seen from the S functions in (D7), in the inter¬ 
esting for us region utx- Ea{ 0) — Ep{0) there is only 
the contribution of the lowest K~ spectrum branch. 

Further evaluation is easily done taking p — 0, for sim¬ 
plicity. Thus, we obtain that for positive e the imaginary 
part differs from zero only for e > lok- (0) + Ap(0) and is 
equal to 


ImE^(e, 0) 


2 TO A 


47r a 


Zq r2(g)g36l(g-pF.p), (D8) 


where g is a solution of the equation A_(e, 0, g) = 0 and 
a = —dA-{e,0,q)/dq ^|q 2 ^g 2 >0. If g is small, namely. 


-A_(e,0,0)^^^M 


q‘^—0 


<C 1 , we may approximate 


g « go = A/A_(e,0,0)/ao, where ao = 

^ q^—0 

As it follows from the 6 function in (D8), fluctuations 
contribute to the imaginary part of the given diagram 
only, if e > UK^iO) + EpiO) -l-p|_pao- 

From Fig. 16 we see that W/f- (q) is very flat function 
of q^, and we can, therefore, neglect the momentum de¬ 
pendence of the K~ spectrum for 0 < jqj 2 TO,r- Then 
we can estimate the imaginary part of the diagram at 
finite p as follows 


Having in mind that we are interested in the self-energy 
insertion to the A — p~^ loop, expanding denominators 
of the fermionic Green’s functions near the poles, and as¬ 
suming that fermions are nonrelativistic, we may put (p ■ 
q) ~ rrfjqu, p^ ~ m*^ and w, jqj, (to^ - to^) < mX,m%. 
Then we come to the non-relativistic p-wave vertex 
Vi ~ i;gq2, with wg = ~ (0.3-b 0.5 )/to2 . 

For E hyperon we would have Vq = '^*n /— 

(0.07 -b 0.09 )/toJ, and for E* the corresponding coupling 
is I’d = I ~ (0.3 -b O.sj/m^. 

To evaluate ImE^ and ReE^ we use the quasiparti¬ 
cle approximation for the spectral functions (B9). As 


„ Ip 1+90 

ImE^(e,p) _ _ 2 f djqj jqp r^- o 

2ml ° J 87t\p\ 

l|p|-9ol 

x2[TO^ + A_(e,O,O)]0[A(e,pF.p,O)] 

« T\qo)qo (g1 +p2)0[A(e,pF.p,O)]. 

(D9) 

The energy and momentum of the A within the Ap~^ 
loop contributing to the K~ polarization operator at 
the zero external momentum q are e = Epip) 4- ux-iO) 
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with IpI < pf,p- Therefore the condition of non-zero 9- 
functions in (D8, D9) is not fulfilled within the momen¬ 
tum integration interval of the Ap~^ loop . 

For ReS-^ using (B5), (B9), (D3), (D6) from (D2) we 
get 


ReY,^{e,p) 

2 ml 


E 


Q 2 2t~i2 


(27r) 


' 7 

'^p—q 


K' 




Ah(e,P,g; 


^q 

A*+(e,p,q 


t}(D10) 


For relevant e the main contribution is given by the first 
term with the lowest K~ branch. For |p| = 0 the integral 
in (DIO) is determined by |q| ~ g, and we can expand 
A_(e, 0, q) « a — q^). If we have 


In order to estimate the Lindhard’s function here, we use 
(33), having 


^pA{t^K- (0)) 0) 

^ 2 mpf Pp ^ Pp 

^pAi^K~i0),0,PF,p) A’ 


A = lok- (0) - Ea{ 0) + Ep{pF,p) ■ 


(D13) 


The value of Ea{ 0) — Ep{pF,p) can be estimated from 
Fig. 6 where it is shown by the dash-dotted line. Us¬ 
ing it, we estimate A ~ —m-Tr and therefore d)pA ~ 
—0.5m^ Pp/po- Thus, we obtain r(0) ~ 1/(1-I-0.3 pp/po) 
and 

Rp yfi 

< 0.2m. (D14) 

2mA 


dZE^ 


dq^ 


q=q 


2 ryK i ^T 

Pf,p ‘E Zg , ^^2 


q=q 


^Pf,p ‘E r(q). 


the remaining integration is straightforward, and we ob¬ 
tain 


ReS^(e,0) _ ^gr"(g)Zf 
2 ml 2 ’k'^ a 


X 


Pf,p 


pf,p 




kl +PF,p 

\q\ - PF,p 


(Dll) 


The extension for finite p can be easily done if we ne¬ 
glect at all the q dependence of the kaon spectrum, that 
is minor, as it is justified by our numerical analysis for 
|q| ;6 2m.. Then we can approximate A_(e,p — q,0) « 
{Qo ~ {P ~ and the integration gives 


ReS"(e,p) 
2 ml 


■r^{qo)z, 


{%+P ) kol 2 In 


^ + (9o +P^)PF.p 


+ PF,p 


ko| -PF,p 


■ (D12) 


Let us estimate the real part of the hyperon self-energy 
at p = 0 and e = Ep{0) -|- lok- (0). We have then qo ^ 0 
and 

» c = Uq mh F^ {0)Zq Pp ~ 0.3m.r^(0) — , 

Z TTlp^ Pq 


for Pp < 3po. We conclude that absolute value of the 
fluctuation contribution is small at p —> 0 and can be 
mimic by a variation of weakly constrained parameters 
of Lagrangian (2). 

From (Dll) and (D12) one can see that ReS^ is log¬ 
arithmically divergent at |p| ^ PF,p when q ^ PF_p. In 
order to analyze, which effects it can induce on the Ap~^ 
contribution to the polarization operator, we separate the 
leading divergent term 


/ReS«\ 

V 2mX 


2™* _ \r 7 K- Pf,p PF,p — IpI 


Vq mh r^(pF,p)2’, 
PF,p - \p\ 


3c In ■ 


PF,p 


PF,p 

(D15) 


and estimate the variation of the Lindhard’s function 
(28): 




PP,P 

f d|p| IpP 

J tt^ 


1 

A - 3c ln(l -p/pF.p) 


0 

Using estimation (D14) we write 


1 

A 



= a 



1 — a In X ’ 


0 


(DI 6 ) 


and find that 


(54>pA 

< 0.2 for 

3c 



A 


where we have used m^ ~ 0.6 mA? and Z^ = ZE^q ~ 
1/(20;^:-(0)) ~ l/(2pe). Assuming that the modi¬ 
fication of the p-wave KNA vertex is determined by 
the graphical equation (71) with the Ap~^ intermediate 
states, we can use the corresponding part of (73) and 
write 

" l-/)vCo<i>pA(u;^-(0),0) ■ 


Thus, the above estimations prove that the self-energy 
corrections of the A propagator induced by the kaon fluc¬ 
tuations (diagram (Dl)) do not modify properties of K~ 
excitations with small energies and momenta, which we 
consider in the main part of this paper. The same esti¬ 
mation can be done also for K^n contribution to the A 
self-energy and for S and S* self-energies. Note that the 
results (D14,D15) and (D16) are the density-independent 
estimations from above. 
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2. Correction of the Baryon Green’s Functions Due 
to Pion Fluctuations 


There is another type of diagrams 
TT 



(D17) 


relating to pion fluctuations. Pions are rather soften in 
nucleon matter already at densities ~ po- This softening 
effect may result in the pion condensation at densities 
p > PcTT, where the critical density for the pion conden¬ 
sation, pcTT, might be smaller than that for the kaon con¬ 
densation [24, 66]. However, concentrating in this work 
on the analysis of the kaon polarization effects we disre¬ 
gard such a possibility assuming p < pcn- 

Most essential contribution to (D17) comes from the 
width of soft pions due to the Landau damping. Max¬ 
imum of the pion spectral function is achieved at small 
pion energy qg < and finite momentum jqj = for 
soft pions. In the diagram with intermediate states 
typical momenta of the hyperons are ~ PFn, if this self¬ 
energy insertion enters the Sn“^, parts of the 

K~ polarization operator, or they are ~ pFp, if this self¬ 
energy insertion enters the Tt*p~^ parts of 

the K~ polarization operator. In a simplifying assump¬ 
tion |q,^| <C 2pFn, 2pFp the nucleon Green’s function can 
be factored out from the integral. The result is then re¬ 
duced to the calculation of the pion tadpole [67] 


vr 


( 

H\ /H 


(DI8) 


At zero temperature the in-medium contribution of 
pion fluctuations given by (DI8) is numerically small 
except for a narrow region near the pion condensation 
critical point, cf. [24, 67]. For finite temperature such 
contributions are substantially increased in the vicinity 
of the pion condensation critical point [68]. 


3. Fluctuation Contributions to the K 
Polarization Operator 


The fermion self-energy insertions discussed above, en¬ 
ter corresponding loop diagrams of the kaon polariza¬ 
tion operator. Since we argued that these contributions 
are rather small we may still work with the quasiparticle 
fermion Green’s functions in the loop diagrams treating 
fermions on the mean-field level, as we did it in main part 
of the text. 

One of important fluctuation processes in the kaon po¬ 


larization operator is given by the diagram 

K 



(D19) 


The diagram with free vertices was under extensive dis¬ 
cussion [19, 20] in application to kaonic atoms. Both the 
vertices were for simplicity taken constants [19], being 

equal to Vo 47r(l -|- mK/mN)^J{a^^y + 3(a^^)2 for 

isospin symmetrical case, where is the KN scat¬ 
tering length for isospin I. In the threshold region one 
obtains a large value of the vertex Vq ~ dTr/m-Tr. We 
will also use the energy independent bare vertex (hatched 
box) but estimate Vq using our amplitudes of Figs. 2, 3 
in the far off-mass shell region. As we see in this energy 
region Vo ~ i.e. a much less value than that in 

the threshold region. Already due to this we expect sig¬ 
nificant suppression of the contribution of this diagram 
to the kaon polarization operator. Another suppression 
comes from short-range correlations taken into account 
in the full box (right vertex). 

With the help of relations (B2) from (D19) we obtain 


m'^(fc) = J ]^-ir+-{k + q)D^{q) 

+V^{k + q) iD^~{q) — V^~(k -I- q)2 Iml?"^(q) 


+iV^{k + q)tD^{q) 


d^q 

(27r) 


4 > 


(D20) 


where V corresponds to the NN~^ loop (including ver¬ 
tices). The last term in (D20) is equal to zero. For fur¬ 
ther convenience we use here the advanced rather than 
the retarded quantities. 

Let us consider first the real part of the diagram. Using 
(B5), (B6), we obtain 


SReU^ik) 


[1 -V n^{ko - qo)] ivik - q)ReD^{-q) 


-f [I -f n^{qo)] AK{q)ReV^{k + q) 


d^q 

(^’ 


(D2I) 


where y-p = 2Im7^'^ and in the first term we have re¬ 
placed q ^ —q. Using relations I -I- = —n^g and 

(D3) we reduce integration in (D21) to the positive en¬ 
ergies: 


(5Ren^(/c) = J 


d^q 


6'(qo) 


[l + n^{ko - qo)hv{k - q)ReD^.{q) 

+ [1 + {ko + qo)]jv{k + q) Re D^+ (q) 

+ [l + n^{qo)] Aji+ {q)Re'P^{k + q) 
+n^{q^)AK-{Q)^eV^{k-q)]- (D22) 
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The NN~^ loop is suppressed at large frequencies uj > 
Therefore, the second and the third 
terms give a small contribution in the energy region of our 
interest. The fourth term is identically zero {n^(qo) = 0 
for <70 > 0 at r = 0). Thus we may keep only the first 
term in (D22) and write 

6ReU^{k) « J ^0{q,)Vo^r^[l + n^{ko-qo)] 

xjo{k-q)^eDi^{q). (D23) 

Here we introduced j-p = T^ 70 using the structure 
of the correlation factor T being expressed via the Lind- 
hard’s function that leads to the T^ factor. In the 
space-like region the 70 (fc) = Joiko,k) has a compact 
analytic presentation, cf. [69], which is further simplified 
at T = 0 of our interest: 


For ko ~ uJx-{0) > Pp j^/{2m*j^) there is always ko > 
and we find 


Re5n^^(A:o,0) 




‘■N 


dTT^ 

C pf,n 
47r^ 


r2(o)4 


PN PF,N 

(D28) 


For ko ~ ujk-{ 0) < pp , inequalities ko > 

H_(q) or ko < H_(q) together with 0 < [qj < 2 pF,n 
determine three regions of the momentum integration. 


Re5nf(fco,0) 

[2 Try 



q+{ko) 


g-(fco) 


2 TT^ 


7o(fco,fc) 


m*j^ ko 
27r|fc| 


, 0 < fco < H_(fc), 


(D24) 


and 


7o(fco,fc) 


* 3 

TYl ^ 

^^iSi_ik) + ko){n^{k) 
0 < H_(fc) < ko < H+(fe), 

^ ^ 2m% 


ko ), 


(D25) 


Let us discuss the low energy contribution (D24) to the 
polarization operator which we label below by subscript 
“1”. Integration in qo is easily done using that 

max(i/, 0) < qo < ko , IT = ko — [k — q). 

Then we get 


ReSU^ik) 


Vq rrij^ 


f d^q T^Zf 

|fc-q|<2pp,„ (27r)3 |fc-q| 

(ko-ujK-iq)) 


(27r)2 
(fco - Wjf-(q))ln 


max(:^, 0) — iVx- (q) 


— ko + max(i/, 0 ) 


(D26) 


This contribution can be easily estimated for \ko — 
'^if-(0)l/^o ^ Wif-(O) and fc = 0. We note that due 
to the flatness of the kaon spectrum we have (q) ~ 
uiR- (0) under the integral. Hence, ko is near the branch 
^K- {q ) we can neglect the first term in (D26). Then we 
have 


Rem7(fco,0)cs-^|^ 

2 Pf,n 

X J min(fco,H_(q)).(D27) 

0 


bordered by q±{ko) = pf,n ± ^Pf^n ~ 2 w^fco satisfying 
ko = fl-{q±). Integration gives 


Re5nf(fco,0) 


C pf,n 

dTT^ 


X 


1- I- 


2 k, 


3/2 


0 Ttlj^ 


Pf,n 


(D29) 


Now we discuss a contribution to the polarization op¬ 
erator from the energy region (D25) which we indicate 
by subscript 2 . 


2pF,iV 




T^Z^^ e{ko-n.{q)) 


0 


.fco-n, (q ) (g ) + fcp _ gp) ^q)_ko + qo) 

>ko-Fl+(q) qo-^K~ 

For ko ~ ujk- (0) and fc = 0 it contributes as 


‘^PF,N 

Re5n2(fco,0)_ ^^4 j |g| 


F^ Z^^ 


x6»(fco - H-(g) ^+{q) In 

-^n^__{q) + 2n^{q)n+{q) - ]^n\{q) 


We have neglected the terms cx {ko — 0JK-{q)) under 
integral, which are small due to the flatness of the K~ 
spectrum. For ko ~ (0) > Pf Ar/(2wJ^) we find 


Re(5n7(fco,0) ~3^^^G(0), (D30) 

and for ko ~ ujx- (0) < Pfn/(‘^''^n) have 
Re SU^ (fco, 0) - 3 dx ), 

x= 2kom*j^/p^j^ , (D3I) 
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where we introduce the function 


duce integration in (D34) to the positive energies: 


G{x) = 



- +4t 


(1 + a:)(3 -x)'^ 1 + 

In ■ 
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3 — a; 


(1 — a:)^(3 + 1 — a; 

— In ■ 


64 


3 + a; ’ 


G(0) = G{x^l)^-^{l-x). (D32) 


The suppression of the nucleon-nucleon-hole loop in the 
scalar-isoscalar channel can be taken into account as 
in Ref. [24] (for recent review see Ref. [70], eq. (4 - 
8)) r = 1/[1 — 2{f + f') Co Ann{(^ = 0, g = pf.at)], 
where Ann is the nucleon-nucleon-hole loop (without 
spin degeneracy factor 2), and Ann(!-^^ = 0, g = py,n) — 
—mArpF,iv/(2 7r^). The Landau-Migdal parameters of 
NN interaction are / — 0 and /' ~ 0.5 0.6 [71]. Thus 

we have F ~ 1/(1 -|- 0.3 {pN/PoY^^)- For pN ^ {3 ^ 5)po 
we have C ~ When lon- (0) > (2to^), we es¬ 

timate the attractive contribution of the diagram (D19) 
as 


ImSn^ = Jeiqo)^^ 

-F g)A;^+(g)[n'’(go) - n‘’(fco + go)] 

+lv {k - q)AK- (g) [1 + n(go) -h n{ko - go)] |d35) 


The second term gives the main contribution. Taking 
into account that at T = 0 there is (1 -I- -I- = 

Q{ko — go)) we have 

ImJn^ = - 

2 

fco ^ 

0 

The K~ spectral function has the form 


Ax-iq) = ‘2TTZg 5{qo-ujK-{q)) 

, _ 7 ^)_ 

bo - "i/C - Ren' 4 (g )]2 + 172 (g)’ 


(D37) 


-ReSn^ < 0.3to^. 

For larger densities, when ujn-{ 0) < pp n/^ this 
contribution is attractive. Its absolute value is addition¬ 
ally suppressed by the ratio tox- (0)?^b/7'F n- 

The imaginary part of the diagram under consideration 
describes the processes when a kaon excitation dissolves 
into multi-particle nucleon-nucleon-hole modes. If the 
K~ energy reached the electron chemical potential in the 
region, where there is the width contribution given by the 
imaginary part of this diagram, the K~ condensation 
would not occur. 

Using (D20) we obtain 

ImSn^ = J ^[tV+-{k + q)lmD^{q) 

+ lmV"^{k + q)iD~^~ (q)] . (D33) 
With the help of relations (B5), (B6) we find 


where we separated two different contributions, the or¬ 
dinary i5-function term from the branch and a possible 
term, which can be obtained only by self-consistent solu¬ 
tion of (D36). 

Let us now consider the quasiparticle contribution to 
the spectral function given by the first term in (D37). 


l{ko,k) 


d^g 


♦ 2 ^2 

t-n *^0 


J (27r)3 2 7r |fe - q I [ 

| fc -( j |< 2 pF , N ' 

[fco - UJX- (g)] 0{ko - UJX- (g)) 

xe{uJx-{q) + fl-(fe - g) - fco) 

-0^[n_(k-q) + ko-cox-{q)] 

X [n+{k -q)-ko+ uJx- (g)] 
x6»(fco - fl_(fc - q) - ujx-iq)) 


ImSU^ = J ^J^^'yp{k + q)Ax{q) 

X [n'^iqo) - n^{ko + go)] . (D34) 


xe{n+{k - q) + ujx-{q) 


fco) 


Using relations 1 -I- n°{—qo) = —n°{qo) and (D3) we re- 


Let us, as before, assume that uJx-{q) is a very flat 



38 


function of q^. Then we obtain 


7 = 


Pf,n \ Vq 


<?+(z) 




rr?*^ T/2 

Stt^ 


('N 
‘^PF,N 


47^3 


^fr^kldlQl 




-q-i-z 


j ^zfr»|n_(<,) + i||n+(,)-2] 


N 

X ° Z r^PF,JV-\/p| AT - 2TO5^f-2: 


27r3 
+0(z)0( 


/I 

4 —i-z 

^^r 
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3^7r3° Pf,n H{2m% z/pp j^), 

2 

Hix)= J ^[4P-iP-xr] 


/T+x-1 


= (! + - + \/l + a:)^ - + a;^ In 


vT+x — 1 


i7(x^0)«4(l + x), i7(x^8)« -(x-8)^(D38) 

o 


here z = ko — ojk- (0). As follows from the 0-functions 


there is no width for energies ko < Ta!aiq{ijJx- (g)} = 
ujR-iQm) for ko < ujk-{^)- The width exists only 
above the branch. Exactly at the branch, i.e. for 
ko = Wat- (fc), the ^-function contributes only for k ^ 0. 
Thus the quasiparticle part of the K~ spectral function 
generates no width at fc = 0 at the critical point of the 
second-order phase transition to the s-wave K~ conden¬ 
sation: ko = Pe = ^K- (0))- 

In the region of low energies, ko < ujx-{k), where 
the (5-function term does not contribute, there might ap¬ 
pear another contribution from the self-consistent solu¬ 
tion given by the second term of (D37). 

Let us find this self-consistent solution and demon¬ 
strate that the width exists even for fc = 0 in this case, 
thus affecting critical condition for the s-wave condensa¬ 
tion. In order to avoid rather cumbersome expressions, in 
our estimate we will do several simplifying assumptions. 

We will assume that ko ~ oJr- (fc) ^ ^FAf, and 

is very small in the interval 0 < |q| < ko/vpN- Then the 
spectral function simplifies as 


^Kiqo) 


_ {Zq )S(go) _ 

[q,-coKA0)?+uzrii<io)y 


(D39) 


and the self-consistent solution of equation (D36) for the 
width 7 is determined by 
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liko) 




2pF,N 


min{/to ,n_ (q)} 


Stt^ 


J d|q||q| J 
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z-f{ko - z) 


(r7K-\'^Y2 f d|g| 

■ 16 7r4 ^ J \q\ 
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0 f2-(q) 

q-(ko) 2pF,iv 
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(fl„(g) -t z) (fl+jq) - z) djko - n^{q))j{ko - z) 
{ko- z- ujk- (0))2 -t \ {Z^~ 7(fco - z)y 


0 q+iko) 


S2-(q) q+{ko) ko 

d|g||g| f dz+ f d|g||g| f dz 

0 q-{ko) 0 


Z7(fco - z) 


{ko-z-ujK-{dW + j{Z^ 'y{ko-z)y 


'>^*N^0 (ryK-\'^T^2 


q-(ko) 2 pF , N ' 


0 q+iko) 


d|g | 

\q\ 


min{feo,S2F(q)} 


dz- 


n_(q) 


(fl_(g) +z) (0+(g) - z)-f{ko - z) 
{ko- z- UK- (0))2 -p i 7(fco - z)y 


(D40) 


where we did the replacement qo = ko — z. For ko ~ uk- 'C ^/2m^ we have g_ « fcoW^/PF.Af) and g_|_ 
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2pF,iv — ^0 iti’^/pf,n ■ Then (D40) reduces to 






‘^Pf,n 


\ iql PF,N _ 

\ 2mj^ 


|g|d|q| 


|g|d|g| 
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kQTrij^ 

2pF N - — 

DTn Rr 


; 2pF,Ar ko 

dz+ J d|q||q| J dz 


z-f{ko - z) 


{ko - z - uiK~{0))^ + J {Z^ -f{ko-z)Y 
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I ^0 ( ryK-\"^ 
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dz- 


(»^(g) + z) (»+(g) - z)7(fco - z) 

{ko- z- ujk- (0))2 + 3 {Z^~ 7(fco - ^)) 


2 ■ 


To solve this equation we assume ^{ko) — ceko for small 
values of ko which we consider here. The main contribu¬ 
tion comes from the third integral of the first line. It is 
(X ko- In the first integral of the first line typical values 
|q| ~ ko and z ^ ko- In the second integral of the first 
line 2pF,Ar— |q| ~ ko and again z ox ko- Then we estimate 
the first two integrals of the first line as ex k^- In the first 
integral of the second line with the help of the replace¬ 
ment z = \q\pF,N/m*jq — q‘^/{2m%) -|-^ we easily see that 
^ ^ q^ ^ k^ and the integral is oc /cq. In the second 
integral of the second line besides this replacement we 
introduce the replacement |q| = 2pf,n — kom’^/pF,N + y 
and observe that ^ cm y (x ko- Then we come to an esti¬ 
mation that the integral is cx fcg. Thus, remaining only 
the third integral of the first line we obtain 


ako — 


'-N >'0 




47r' 


' 2 r2 „2 

^ Pf,n 


4 

za{ko — z)dz 


{ko - z - ujk-{ 0W + j{Z^ ) a'^{ko-z)'^ 


which has a non-trivial solution: 


7 ~ ako - rnyVoTpF^Nko/ir'^ , (D4I) 


for a > 2/Z(^ . As one may expect, this inequal¬ 
ity indeed holds for rather large densities since Z^ ~ 
Xj[2iOx-{^)) for small values ko ~ ujx~{0)- 

The possible presence of a width at low energies is 
a principal question. The second-order phase transi¬ 
tion to the s-wave K~ condensation cannot occur, if 
tOfc- (0) crosses pe in the energy region, where exists 
imaginary part of the polarization operator. The pos¬ 
sibility of the condensation is however recovered, when 


the electron chemical potential exceeds lok- (0) reach¬ 
ing the upper boarder of the region of the width. At 
least for lok- (0) > Pfn/ not find self- 
consistent solution for the width. Thereby we would 
like to stress that the second order phase transition 
to the s-wave K~ condensation may not occur only if 
zJK{d) = pe < p|„/(2m^). For realistic values of pa¬ 
rameters we have (0) = pe > Ppm/ (2w^) without any 
problem for the s- wave K~ condensation. 

Thus, by the above examples we demonstrated that in 
spite of many new peculiarities associated with the fluc¬ 
tuation processes for a rough analysis we may drop the 
fluctuation contributions treating baryons on the mean 
field level. 

Similar diagram to (D19) (but with the pion instead of 
the kaon and a different vertex) also describes the K~tt — 
ttK~ interaction. In a simplifying assumption that the 
pion momentum |< 7 „| <C 2pF,n the diagram is reduced to 

vr 

(D42) 



the same as (D18) but with a different vertex and the 
kaon external Green’s functions. As we mentioned such 
a contribution can be easily estimated. At zero tem¬ 
perature the in-medium contribution of pion fluctuations 
given by (DI8) is numerically small except for a nar¬ 
row region near the pion condensation critical point, cf. 
[24, 67]. For finite temperature such contributions are 
substantially increased in the vicinity of the pion con¬ 
densation critical point [68]. 
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